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(54) Title: SYSTEMS AND METHODS FOR CONSTRUCTING GENOMIC-BASED PHENOTYPIC MODELS 



(57) Abstract: The invention provides a computer implemented process for constructing a scalable output network model for a 
bioparticle. The process includes computer implemented steps of: (a) accessing a database of network gene components includ- 
ing an annotated network set of open reading frames (ORFs) of a bioparticle genome, (b) forming a data structure associating the 
network gene components with network reaction components, the data structure establishing a data set specifying a network model 
of connectivity and flow of the network reaction components, and (c) transforming the data set into a mathematical description of 
reactant fluxes defining the network model of connectivity and flow, wherein the mathematical description defines a scalable output 
network model of a bioparticle. 
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SYSTEMS AND METHODS FOR CONSTRUCTING GENOMIC -BASED 

PHENOTYPIC MODELS 

BACKGROUND OF THE INVENTION 

This invention relates generally to simulation 
5 modeling and, more specif ically, to computational methods 
for simulating and predicting the activity of biochemical 
and biological network models. 

Therapeutic agents , including drugs and 
gene-based agents, are being rapidly developed by the 

10 pharmaceutical industry with the goal of preventing or 
treating human disease. Dietary supplements, including 
herbal products, vitamins and amino acids, are also being 
developed and marketed by the nutraceutical industry. 
Because of the complexity of biochemical reaction 

15 networks, even relatively minor perturbations caused by a 
therapeutic agent or a dietary component on the abundance 
or activity of a particular target, such as a metabolite, 
gene or.. protein, can affect hundreds of biochemical 
reactions. These perturbations can lead to desirable 

20 therapeutic effects, such as cell stasis or cell death in 
the case of cancer cells or other pathologically 
hyperprolif erative cells. However, these perturbations 
can also lead to undesirable side effects, such as 
production of toxic byproducts. 

25 Traditionally the identification of drugs and 

nutraceuticals has relied upon early stage screening and 
testing in which the effects of candidate drugs on 
individual genes or gene products are observed. This 
approach, although helpful for identifying a particular 
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gene or gene product as a target for a particular 
disease, is often incapable of identifying the effects 
that the candidate drug or the drug inhibited target will 
have on other molecular components of the cell or 
5 organism. It is often not until late stage testing with 
human subjects that unwanted or even dangerous side 
effects are observed. Failure to select a candidate drug 
in early stage testing that is without side effects can 
result in harm to individuals participating in clinical 
10 trials and significant delays in curing individuals 

suffering from disease due to pursuing the wrong drug. 

In order to design effective methods of 
repairing, engineering or disabling cellular activities, 
it is essential to understand cellular behavior from an 

15 integrated perspective. Methods have recently been 

developed to reconstruct biological reaction networks 
that occur within organisms, with the goal of being able 
to model them and then use simulation to predict and 
analyze organismal behavior. One of the most powerful 

20 current approaches to modeling complex biological 

reaction networks involves constraints-based modeling. 
This approach provides a mathematically defined solution 
space wherein all possible behaviors of the reconstructed 
biological reaction network must lie. The solution space 

25 can then be explored to determine the range of 

capabilities and preferred behavior of the biological 
system under various conditions. 

A combination of many high throughput 
technologies is now providing information on a scale that 
30 includes entire genomes, the complete set of gene 

products encoded by the genomes, and molecular functions 
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that occur in a cell or organism. The ability to create 
genome scale constraints-based models requires that vast 
amounts of biological information be assimilated. 
Although genome scale models have been produced for a 
5 variety of organisms and have been shown to accurately 
predict a number of cell functions, it is currently 
difficult and time consuming to build new models and many 
organisms for which genome scale information is available 
currently lack genome scale models. Furthermore, it is 
10 currently difficult to view the content of models and to 
cross-reference the information in the models with the 
information available in biological databases and with 
other models. Thus, for many models, errors go unnoticed 
or are difficult to correct once the model is built. 

15 

Thus, there exists a need for constraints-based 
models for the increasing number and variety of organisms 
for which genomes are being sequenced. A need also 
exists for methods to efficiently build and modify 
20 existing constraints-based models. The present invention 
satisfies these needs and provides related advantages as 
well . 

SUMMARY OF THE INVENTION 

The invention provides a computer implemented 
25 process for constructing a scalable output network model 
of a bioparticle. The process includes computer 
implemented steps of: (a) accessing a database of network 
gene components comprising an annotated network set of 
open reading frames (ORFs)of a bioparticle genome; (b) 
30 forming a data structure associating the network gene 
components with network reaction components, the data 
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structure establishing a data set specifying a network 
model of connectivity and flow of the network reaction 
components , and (c) transforming the data set into a 
mathematical description of reactant fluxes defining the 
5 network model of connectivity and flow, wherein the 
mathematical description defines a scalable output 
network model of a bioparticle. 

The invention .further provides a computer 
implemented process for constructing a scalable 
phenotypic output network model. The process includes 
the computer implemented steps of: (a) accessing a 
database of network gene components including an 
annotated network set of open reading frames (ORFs)of a 
bioparticle genome; (b) forming a data structure 
associating the network gene components with network 
reaction components, the data structure establishing a 
data set specifying a network model of connectivity and 
flow of the network reaction components; (c) modifying 
the data set to enumerate a biochemical demand on the 
specified network model/ and (d) transforming the 
modified data set into a mathematical description of 
reactant fluxes defining the network model of 
connectivity and flow, wherein the enumerated biochemical 
demand corresponds to an aggregate reactant demand flux 
defining a phenotypic output of the network model of a 
bioparticle . . 

Also provided is a computer implemented process 
for self -optimizing a network model of a bioparticle. 
The process includes the computer implemented steps of: 
30 (a) accessing a database of network gene components 

including an annotated network set of open reading frames 
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(ORFs)of a bioparticle genome; (b) forming a data 
structure associating the network gene components with 
network reaction components, the data structure 
establishing a data set specifying a network model of 
5 connectivity and flow of the network reaction components; 
(c) transforming the data set into a mathematical 
description of reactant fluxes defining the network model 
of connectivity and flow; (d) determining the competence 
of the connectivity and flow within the network model, 

10 the competence indicating underinclusion or overinclusion 
of network reaction component content of the network 
model, and (e) identifying an ameliorating network 
reaction component capable of augmenting the competence 
of the. network model, incorporation of the ameliorating 

15 network reaction component into the data structure 

producing a modified data structure specifying in an 
optimized network model of the bioparticle. 

The invention also provides a computer 
implemented process for constructing a data structure 

20 specifying a network model of a bioparticle. The process 
includes the computer implemented steps: (a) accessing a 
database of network gene components including an 
annotated network set of open reading frames (ORFs)of a 
bioparticle genome; (b) selecting an ORF from the 

25 annotated network set encoding a gene product having a 

network reaction function; (c) determining the occurrence 
of a constituent gene product for the selected encoded 
gene product; (d) determining the occurrence of an 
additional gene product participating in the network 

30 reaction; (e) forming a data structure from the selected 
and determined gene products, the data structure 
associating the network gene components and network 
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reaction components comprising cognate ORFs, encoded gene 
products,, network reactions and reaction constituents, 
and (f)' repeating steps (a) -(e) selecting another ORF 
from the annotated network set until substantially all of 
5 the network gene components of the annotated network set 
have been surveyed for encoding a gene product having a 
network reaction function to produce a data structure 
establishing a data set specifying a network model of 
connectivity and flow. The invention further provides 
10 computer systems having executable instructions for 
carrying out these computer implemented processes. 

A system for constructing a scalable output 
network model of a bioparticle, including: (a) an input 
data set of network gene components including an 

15 annotated network set of open reading frames (ORFs) of a 
bioparticle genome; (b) executable instructions forming a 
data structure- associating the network gene components 
with network reaction components, the data structure 
establishing a data set specifying a network model of 

20 connectivity and flow of the network reaction components; 
(c) executable instructions determining the occurrence of 
a reaction component satisfying a macro requirement 
deficiency in structural architecture of the network 
model, inclusion of an identified reaction component 

25 satisfying the macro requirement deficiency in the data 
structure supplementing the connectivity and flow of the 
network model; (d) a heuristic logic decision algorithm 
determining confidence of the network reaction components 
within the data structure, and (e) executable 

30 instructions mathematically describing from the data set 
reactant fluxes defining the network model of 
connectivity and flow, wherein the mathematical 
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description defines a scalable output network model of a 
bioparticle. A system for constructing a scalable 
phenotypic output network model of a bioparticle, 
including: (a) an input data set of network gene 
5 . components including an annotated network set of open 
reading frames (ORFs) of a bioparticle genome; (b) 
executable instructions forming a data structure 
associating the network gene components with network 
reaction components, the data structure establishing a 

10 data set specifying a network model of connectivity and 
flow of the network reaction components; (c) executable 
instructions modifying the data set to enumerate a 
biochemical demand on the specified network model, and 
(d) executable instructions mathematically describing 

15 from the modified data set reactant fluxes defining the 
network model of connectivity and flow, wherein the 
enumerated biochemical demand corresponds to an aggregate 
reactant demand flux defining a phenotypic output of the 
network model of said bioparticle. A system for 

20 constructing a self -optimizing network model of a 

bioparticle, including: an input data set of network gene 
components including an annotated network set of open 
reading frames (ORFs) of a bioparticle genome; executable 
instructions forming a data structure associating said 

2 5 network gene components with network reaction components, 
said data structure establishing a data set specifying a 
network model of connectivity and flow of said network 
reaction components; executable instructions 
mathematically describing from said data set reactant 

30 fluxes defining said network model of connectivity and 
flow; executable instructions computing competence of 
said connectivity and flow within said network model, 
said competence indicating underinclusion or 
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overinclusion of network, reaction component content of 
said network model, and executable instructions 
augmenting said competence of said connectivity and flow 
within said network model, said executable instructions 
5 specifying inclusion or exclusion of an ameliorating 

network reaction component, wherein incorporation of said 
ameliorating network reaction component into said data 
structure produces a modified data structure specifying 
an optimized network model of said bioparticle. 

10 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 shows an exemplary system architecture 
for a computer system of the invention. 

Figure 2 shows an overview of an exemplary 
model construction process. 

15 Figure 3 shows an associated object model of a 

network model specifying the participating classes of 
network component data elements and associations in a 
biochemical network of a bioparticle. 

Figure 4 shows an associated database schema of 
20 a network model specifying the participating tables of 
network component data elements and associations in a 
biochemical network of a bioparticle. 

Figure 5 shows an exemplary process of 
constructing a data structure of network reaction 
25 components. 
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Figure 6 shows an exemplary process of 
evaluating a gene index and creating reaction 
associations . 

Figure 7 shows association diagrams displaying 
5 ORF-protein-reaction associations. 

Figure 8 shows the Model Construction main 
window for a system of the invention. 

Figure 9 shows a model construction window with 
a display of a gene index for a bioparticle. 

10 Figure 10 shows a model construction window 

with in which the AceEF protein is entered into the 
XN Protein" entry field, thereby being associated to the 
b0114 and b0115 ORFs . 

Figure 11 stiows a model construction window in 
15 which gene-protein associations for the AceEF protein are 
displayed visually in a graphical association viewer and 
the requirement for two ORFs to encode the protein is 
represented by an "AND" association. 

Figure 12 shows a model construction window in 
20 which gene-protein-reaction associations for the 

TRANS (pi) reaction are displayed visually in a graphical 
association viewer and the requirement for two ORFs to 
encode the protein is represented by an "AND" 
association . 

25 Figure 13 shows a model construction window in 

which gene-protein-reaction associations for the PYRDH 
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reaction are displayed and different isozymes that 
catalyze the reactions are represented by drawing 
multiple lines between the ORFs and the protein. 

Figure 14 shows a model construction window in 
5 which a protein that is associated with a model is 
displayed in a table. 

Figure 15 shows a model construction window in 
which a protein that is associated with a model and 
displayed in a table is selected for inclusion in a 
10 model. 

Figure 16 shows a model construction window in 
which ORF-protein~reaction associations are visually 
displayed in a graphical association viewer. 

Figure 17 shows a model construction window in 
15 which a protein-reaction "AND" association is displayed 
in a graphical viewer. 

Figure 18 shows a model construction window in 
which a protein-reaction "OR" association is displayed in 
a graphical viewer. 

2 0 DETAILED DESCRIPTION OF THE INVENTION 

Computer systems and computer implemented 
processes for constructing and using a network model of a 
bioparticle are described. In the following description, 
for the purposes of explanation, specific details are set 
25 forth in order to provide a thorough understanding of the 
present invention. Those skilled in the art will 
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understand that the present invention can be practiced 
without these specific details and can be applied to any 
of a variety of related systems. For example, although 
the methods are described in the context of metabolic 
5 reactions it is understood that similar models can be 
made and used for simulation of other network systems 
such as biological regulatory systems, biological signal 
transduction systems and non-biological reaction systems. 

10 In one embodiment, a network model of the 

invention can be used in sllico to simulate the flux of 
mass, energy or charge through the chemical reactions of 
a biological system to define a solution space that 
contains any and all possible functionalities of the 

15 chemical reactions in the system, thereby determining a 
range of allowed activities for the biological system. 
Such an approach is referred to as constraints-based 
modeling because the solution space is defined by 
constraints such as the known stoichiometry of the 

20 included reactions as well as reaction thermodynamic and 
capacity constraints associated with maximum fluxes 
through reactions. Using a network model of the 
invention, the space defined by these constraints can be 
interrogated to determine the phenotypic capabilities and 

25 behavior, of the biological system or of its biochemical 
components. Analysis methods such as convex analysis, 
linear programming and the calculation of extreme 
pathways as described, for example, in Schilling et al . , 
J. Theor. Biol. 203:229-248 (2000); Schilling et al., 

30 Biotech. Bioena. 71:28*6-306 (2000) and Schilling et al . , 
Biotech, Proa. 15:288-295 (1999), can be used to 
determine such phenotypic capabilities. 
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In another embodiment, the constraints-based 
method is flux balance analysis. Flux balance analysis 
■is based on flux balancing in a steady state condition 
and can be performed as described in Varma and Palsson, 
5 Biotech. Bioencr. 12:994-998 (1994). Flux balance 
approaches can be applied to .reaction networks to 
simulate or predict systemic properties of adipocyte 
metabolism as described in Fell and Small, J. Biochem. 
138:781-786 (1986), acetate secretion from E . coli under 

10 ATP maximization conditions as described in Majewski and 
Domach, Biotech. Bioencr . 35:732-738 (1990) or ethanol 
secretion by yeast as described in Vanrolleghem et al., 
Biotech. Prog. 12:434-448 (1996). Additionally, this 
approach can be used to predict or simulate the growth of 

15 E . coli on a variety of single-carbon sources as well as 
the metabolism of H. influenzae as described in Edwards 
and Palsson, Proc, Natl. Acad, Sci. 97:5528-5533 (2000), 
Edwards and Palsson, J. Bio. Chem. 274:17410-17416 (1999) 
and Edwards et al., Nature Biotech. 19:125-130 (2001). 

20 Once the solution space has been defined, it 

can be analyzed to determine possible solutions under 
various conditions. This is an approach that is 
consistent with biological realities. Biological systems 
have built in flexibility and can, therefore, reach the' 

25 same result in many different ways. These systems are 
designed through evolutionary mechanisms that have been 
restricted by fundamental constraints that all living 
systems must face. The constraints-based modeling 
strategy embraces these general realities. 

30 For a reaction network that is defined for a 

particular organism through the use of genome sequence 
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and biochemical and physiological data, the solution 
space describes the functional capabilities of the 
organism as described, for example, in WO 00/46405. 
Genome scale models have been created for a number of 
5 organisms including Escherichia coli (Edwards et al., 
Proc. Natl. Acad. Sci. USA 97:5528-5533 (2000)), 
Haemophilus influenzae (Edwards et al., J. Biol. Chem. 
274: 17410-17416 (1999)), Bacillus subtilis and 
Helicobacter pylori. 

10 The ability to continuously impose further 

restrictions on a network model via the tightening of 
constraints results in a reduction in the size of the 
solution space, thereby enhancing the precision with 
which physiological performance or phenotype can be 

15 predicted. This approach provides a basis for 

understanding and ultimately predicting the structure and 
function of a biological system through the model 
building and implementation process as set forth below. 

As used -herein, the term ^scalable" is-" intended' 
20 to mean that the content size of a network model . of the 
invention can increase without substantial diminution in 
model performance where performance is a measure of model 
predictability. In general, the performance of a network 
model will increase proportionally to the accuracy of 
25 content elements included in the model. Although the 
number of calculations can increase with increase in 
content size, the predictability for obtaining a 
particular solution for a scalable network model of the 
invention will not be substantially diminished due to 
30 changes in content size alone. Network model content 
that can be increased includes, for example, data 
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elements specifying gene component and network reaction 
components." The scalable network models of the invention 
also includes, for example, increasing network model 
content from a simple system of gene and network reaction 
5 components to complex, multisystem gene and network 
reaction components, to network gene and reaction 
components specifying complex cell and multicellular 
systems without substantial diminution in model 
performance. A specific example of maintaining network 

10 model performance while increasing model content would be 
increasing the model content of a gene to that specifying 
substantially all biochemical reactions derived from a 
cellular genome. Therefore, the term includes the 
ability of a network model to expand the number of ORFs, 

15 reactions, reactants and fluxes without requiring 
manipulations to the model programming, design or 
software architecture . 

As used herein, the term "bioparticle" is 
intended to mean a biological entity that contains a 

2 0 nucleic acid genome that. encodes constituent parts of the 
entity. The nucleic acid genome can be, for example, DNA 
or RNA. and can be derived from a naturally occurring 
biological entity, a non-naturally occurring biological 
entity or designed de novo. A biological entity included 

25 in the term can be, for example, a virus or a cell, such 
as a procaryotic cell or eucaryotic cell or other 
naturally occurring or non-naturally occurring biological 
entities. A cell can be derived from a unicellular 
organism or from a multicellular organism. 

30 As used herein, the term "phenotype," when used 

in reference to a network model, is intended to mean the 
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detectable characteristics resulting from the interaction 
of a model genotype and a model environment. A 
detectible characteristic refers to a computed individual 
or integrated function of one or more network model 
5 components. Network models of the invention simulate, in 
silico, an organism or a functional set of interactive 
components of an organism. A model genotype contains 
those network gene components included in a network model 
specifying an in silico organism. A model environment 

10 includes, for example, a specified external condition 

exposed to an in silico organism. Therefore, a phenotype 
of a network model is a detectable result of the 
functional interactions of gene products encoded in the 
model genotype, and related reaction components, and the 

15 environmental conditions which influence the activity and 
interactions of network model components. A "phenotypic 
output'' as it is used herein, refers to the m'easure of a 
characteristic resulting from simulation of a network 
model, or from simulation of a particular solution to a 

20 network model. A phenotypic output can be, for example, 
a solution space of a network model where the model 
environment consists all possibilities, a feasible 
solution where the model environment consists of 
constrained fluxes of external components, or a 

25 particular solution where the model environment consists 
of defined components. 

As used herein, the term "network" is intended 
to mean a system of interconnected or interrelated 
components. The interconnections and interrelations can 
30 be, for example, either physical or functional 

relationships of system components. Therefore, the term 
refers to an aggregation or assemblage of system 
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components and the relative relationships that define 
inclusion of components within such a system. One 
example of a network can be a computational 
representation of genes, gene products , reactants, 
5 functions and physicochemical characteristics , for 

example, that constitute an in silico organism of the 
invention. Another example of a network can be a 
computational representation of a genes, gene products, 
reactants, functions and physicochemical characteristics, 

10 for example, that constitute a biochemical network or a 
biochemical pathway of an in silico organism. Such 
biochemical networks can include, for example, central 
metabolism, peripheral metabolism, protein biosynthesis, 
carbohydrate biosynthesis, lipid biosynthesis and signal 

15 transduction. Biochemical pathways can include, for 

example, glycolysis, the citric acid (TCA) cycle, amino 
acid biosynthesis, nucleoside and nucleotide 
biosynthesis, a signal transduction event, and the like. 
Numerous, other examples of reactions or events that 

2 0 combine into networks and pathways to produce a common 
function are well known to those skilled in the art and 
are included within the meaning of the term. Such 
networks and pathways can be found described in, for 
example, Stryer, L., Biochemistry , W.H. Freeman and 

25 Company, New York, 4th Edition (1995); Alberts et al., 

Molecular Biology of The Cell , Garland Publishing, Inc., 
New York, 2nd Edition (1989); Kuby, Immunology , 3rd 
Edition, W.H. Freeman & Co., New York (1997), Kornberg 
and Baker, DNA Replication , W.H. Freeman and Company, New 

30 York, 2nd Edition (1992), all of which are incorporated 
herein by reference. Therefore, regardless of the label 
used or the number of constituent elements, a network 
refers to a collection of components that exhibit a 
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logical physical or functional relationship whose 
concerted interaction are employed for at least one 
common purpose. 

As used herein, the term "component" or 
5 "network component" is intended to mean a data element, 
data set or electronic representation of a chemical or 
biochemical molecular entity in a network model of the 
invention. The term is intended to refer to the input 
and output representations as well as to the code and 

10 electronic representations within a computer program or 

processor. Therefore, representations of components of a 
system and their interrelationships will depict a network 
model of the invention. A variety of formats well known 
to those skilled in the art can be used to represent any 

15 or all types of chemical and biochemical components 
within a network model. The term can include, for 
example, a gene component, a reaction component or a non- 
gene component . 

As used herein, the term "gene component" is 
20 intended to mean a data element, data set or electronic 
representation of a nucleic acid that encodes a gene 
product, or functional fragment thereof. A gene 
component can be represented ' in a network model by, for 
example, nucleotide sequence, nucleic acid structure, 
25 name, symbol, . with reference to its encoded gene product, 
activity or combination thereof. The term is intended to 
refer to input and output representations, such as text 
and visual graphics, as well as to programming code or 
electronic representations within a computer processor. 
30 Therefore, a "network gene component" as used herein, 
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refers to a gene component which is part of a network 
model of the invention. 

As used herein, the term "reaction component" 
is intended to mean a data element, -data set or 
5 electronic representation of a component of a network, or 
functional fragment thereof. A network reaction 
component can be, for example, a gene product, a 
macromolecule or a molecule. Specific examples of 
network reaction components include enzymes, substrates, 

10 products, cofactors, DNA, RNA, polypeptide, lipid, 
carbohydrate, amino acids, nucleotides, nucleotide 
triphosphates, fatty acids, sugars, steroids, 
metabolites, catabolites, ions, metals, and the like. 
Such gene products participate or function in a wide 

15 variety of chemical or biochemical reactions well known 
to those skilled in the art, including for example, 
chemical reactions, binding reactions and signal 
transduction reactions. A reaction component can be 
represented in a network model by, for example, primary 

20 structure such as- amino, acid or other monomer sequence of 
a polymer, secondary structure, tertiary structure, name, 
symbol, with reference to its encoding gene, reactants, 
activity or combination thereof. The term is intended to 
refer to input and output representations, such as text 

25 and visual graphics, as well as to code or electronic 

representations within a computer processor. Therefore, 
a "network reaction component" as used herein, refers to 
a reaction component which is part of a network model of 
the invention. 

30 As used herein, the term "network set" when 

used in reference to network gene components is intended 
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to mean a group of network gene components encoding gene 
products that complete a concerted function of a network. 
Therefore, a network set is at least a subset of 
components that constitute a network model of the 
5 invention. A network set also can contain all components 
constituting a network model of the invention. So long 
as a set of components can complete a concerted function 
of a network, a network set can include, for example, 
biochemical networks, biochemical pathways and other 

10 biochemical systems well known in to those skilled in the 
art. A network set is ''annotated" when it is derived 
from a gene sequence record that specifies a function or 
attribute of the recorded gene or a gene product encoded 
therefrom. Because gene records will have at least one 

15 function or attribute associated with them, essentially 
all gene sequences that have been recorded in a tangible 
medium or archived are included within the meaning of the 
term annotated. A function can include, for example, an 
activity of an encoded gene product such as the 

20 conversion of substrate to product or the transition from 
an inactive state to an active state in the presence of a 
stimulus. An attribute can be, for example, a nucleotide 
sequence, a name, a nucleotide or amino acid composition, 
a molecular weight, a size or a structure. Specific 

25 examples of annotated network sets include a genome as 
well as those biochemical networks and biochemical 
pathways exemplified previously with reference to 
networks of the invention. Sources of annotated network 
sets include, for example, Genbank; Unigene; Subtilist 

30 (Bacillus subtilis); YPD (Saccharomyces cerevisiae) ; 
Wormbase (Caenorhabditis elegans) ; ensembl (Human, 
mouse); PKR (kinases); GPCRDB (G-proteins) ; EcoCyc, KEGG, 
WIT, B REN DA (metabolism) ; Regulon DB, Transfac 
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(regulation); and AFC S, TRANS PATH (signal transduction). 
These and other databases from which annotated network 
sets can be obtained are well known in the art as 
described, for example, in Baxevanis, Nucleic Acids Res. 
5 30:1-12 (2002). 

As used herein, the term "data structure" is 
intended to mean an organization of information, such as 
a physical or logical relationship among data elements, 
designed to support specific data manipulation functions, 

10 such as an algorithm. The term can include, for example, 
a list or other collection type of data elements that can 
be added, subtracted, combined or otherwise manipulated. 
Exemplarily, types of data structures include a list, 
linked-list, doubly linked-list, table, matrix, queue, 

15 stack, heap, dictionary and tree. Such organizational 
structures can include, for example, data elements 
representing all categories and subcategories of network 
components. The term also can include organizational 
structures of information that relate or correlate, for 

20 example, data elements from a plurality of data 

structures or other forms of data management structures. 
A specific example of information organized by a data 
structure of the invention is the association of a 
plurality of reactions with corresponding reactants and 

25 stoichiometry for a network model. Other information 
that can be organized by a data structure of the 
invention includes, for example, a representation or 
relationship of a substrate or product of a chemical 
reaction, a chemical reaction relating one or more 

30 substrates to one or more products, a constraint placed 
on a reaction, or a stoichiometric coefficient. 
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As used herein, the term "data set" is intended 
to mean a collection of data elements. A specific 
example of a data set is a file. Hierarchical forms and 
organizations of data sets are also included within the 
5 meaning of the term. Data element refers to a unit of 
data or a computational representations thereof. 
Generally, data elements and data sets are processed or 
interpreted to take on meaning. Data representations can 
include, for example, numbers, characters, images, or 
10 other method of recording well known in the art, in a 
form that can be input into a computer, stored and 
processed there, or transmitted on some digital channel. 
Therefore, data elements can be represented, for example,, 
in machine language, assembly language or user language. 

15 As used herein, the term "connectivity" is 

intended to mean the pattern, interactions and routes of 
linkage between network components. Such linkages serve 
to place network components in a physical or functional 
relationship that specifies a unity of common plan or 

20 purpose of such components. Therefore, the term 

connectivity refers to the aggregation and assemblage of 
network components joined through physical or functional 
interaction or interdependence. For example, a chemical 
reaction that converts compound A to compound B links 

25 these compounds by physical interconversion function 

within a network model. Similarly, where an enzyme uses 
compound B as a substrate to produce product P, the 
enzyme and its chemical reaction is functionally linked 
by interdependence to the above chemical reaction that 

30 produces compound B. A specific example of a complex 

system of connectivity constitutes some or substantially 
all of the biochemical reactions, interactions and 
interdependencies of a bioparticle. 
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As used herein, the term "flux" or "reactant 
flux" is intended to refer to the flow, transfer or 
conversion of a network component through a reaction or 
network. A reaction included in the term can be any 
5 conversion that consumes a substrate or forms a product 
including, for example, changes in chemical composition 
such as those that occur due to an enzymatic process, 
changes in location such as those that occur due to a 
transport reaction that moves a reactant from one 

10 cellular compartment to another or a binding reaction. 
The term includes directionality and can be represented 
by a variety of means and formats known to those skilled 
in the art. For example, conversion of substrate to 
product can be represented as a positive flux of product, 

15 corresponding to its formation; or as a negative flux of 
substrate, corresponding to its disappearance. Positive 
fluxes also can be characterized to have a forward 
direction whereas negative fluxes can be characterized as 
a backward direction. Fluxes also can be represented by, 

20 for example, a reaction showing directionality. The term 
"flux" when used in reference to a pathway or flux 
pathway is intended to include combinations and 
permutations of individual fluxes, such as the flow or 
transfer of network components through a series multiple 

25 reactions. Exemplarily combinations and permutations of 
individual fluxes include a flow, transfer or conversion 
of network components in or through a biochemical pathway 
or a biochemical network. Descriptions or 

representations of a flux or a flux pathway can be either 
30 qualitative or quantitative.. 
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As used herein, the term "aggregate reactant 
flux'" or "aggregate reactant demand flux" is intended to 
mean the combined flow, transfer or conversion of network 
components through reactions of two or more reaction 
5- pathways into a single category for model representation 
or analysis. Combination of reaction pathways can occur, 
for example, at the terminal output of a reaction pathway 
or at any point along the pathway or transfer of 
reactantS: or products. Therefore, an aggregate flux can 

10 be a portion or subset of a reaction pathway. Aggregate 
fluxes can be used to define a variety of external inputs 
and outputs to a system as well as to define internal 
inputs and outputs that are secondary to the primary 
network of a particular model. Therefore, the term also 

15 is intended to include both internal system fluxes and 

external fluxes. For example, an internal aggregate flux 
can be a representation of all amino acid biosynthesis as 
a single reaction flux. An external aggregate flux can 
be, for example, a representation of the import into the 

20 system of all carbon sources used or by-products 

generated in an in sxlico network model of the invention. 
Aggregate fluxes also can be implemented in a network 
model to define the activity of one or more biochemical 
demands . 

25 As used herein, the term "biochemical demand" 

is intended to mean a flux, a flux pathway or an 
aggregate flux that represents a biochemical requirement. 
Such requirements can include, for example, network 
components used for growth or other cellular or 

30 physiological processes, metabolism, catabolism, energy 
production, redox equivalent production, biomass 
production, development, or consumption of carbon 
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nitrogen, sulfur, phosphate, hydrogen or oxygen. 
Examples of a particular network components used for such 
requirements include, for example, the production of 
biomass precursors, production of a protein, production 
5 of an amino acid, production of a purine, production of a 
pyrimidine, production of a lipid, production of a fatty 
acid, production of a cofactor, production of a cell wall 
component or transport of a metabolite. Other 
biochemical demands and their corresponding network 

10 components well known to those skilled in the art also 
included within the meaning of the term- 
As used herein, the term "macro requirement 
deficiency" is intended to mean the absence of flux or 
inappropriate flux directionality from one component of a 

15 network model to another interrelated network component. 
Absence of flux includes, for example, an undesirable 
buildup of a reaction product, lack of a substrate 
required for a reaction to occur, or a gap in a reaction 
network wherein a metabolite can be produced, but not 

2 0 consumed or where a metabolite, can be consumed but not 
produced. Absence or- inappropriate flux also can 
include, for example, singleton network components that 
exist in the system model in isolation and multiple, 
adjacent network components that have irreversible 

25 thermodynamic assignments. • A specific example of a 
singleton network component is a reaction within a 
biochemical pathway existing in a network model without a 
flux of reactants to and from the reaction. A specific 
example of multiple, adjacent irreversible components is 

30 where two or more connected reactions have irreversible 
kinetic parameters . 
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As used herein,, the term "elemental balancing" 
refers to conservation of chemical elements during 
chemical transformation of one network component into 
another. The term therefore includes the stoichiometry 
5 of a chemical reaction as well as accounting for other 
chemical inputs and outputs of a chemical reaction. A 
specific example' of elemental balancing includes ensuring 
that the total number of * oxygen atoms, for example, in 
all reactants used in a transformation equals the number 

10 of oxygen atoms in all the reactants formed by the 
transformation. Similarly, for all other atoms 
constituting the substrates or input reactants in a 
transformation, the number of each type of atom consumed 
will equal the number of the same type of atom formed if 

15 that reaction is elementally balanced. In the case of 
multiple transformations, such as those constituting a 
reaction network, the multiple transformations will be 
elementally balanced when, for each atom, the net number 
of the same type of atom consumed by the multiple 

20 transformations, taken as a whole, is equal to the net 
number of the same type of atom formed by the multiple 
transformations, taken as a whole. Elemental balancing 
includes, for example, all elements within the Periodic 
Table such as carbon, hydrogen, phosphorus, nitrogen, 

25 zinc, magnesium and the like. The term "charge 

balancing" refers to the similar process of accounting 
for equivalent input and output of all electrical charges 
on a reactant participating in one or more chemical 
reactions . 

30 The invention provides a computer implemented 

process for constructing a scalable output network model 
of a bioparticle. The process includes the computer 
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implemented steps of: (a) accessing a database of network 
gene components including an annotated network set of 
open reading frames (ORFs)of a bioparticle genome; (b) 
forming a data structure associating the network gene 
5 components with network reaction components, the data 
structure establishing a data set specifying a network 
model of connectivity and flow of the network reaction 
components, and (c) transforming the data set into a 
mathematical description of reactant fluxes defining the 
10 network model of connectivity and flow, wherein the 
mathematical description defines a scalable output 
network model of a bioparticle . 

A computer implemented process of the invention 
can be carried out on a computer system that provides a 

15 means to construct, access, modify or utilize a network 
model of the invention as well as the information 
associated with the network model. A computer system can 
have any of a variety of known architectures including, 
for example, single tier or multi-tier architectures. An 

2 0 exemplary architecture for a computer system of the 

invention is the multi-tier or multi-server application 
shown in Figure 1 and consisting of an application server 
1 that communicates with a client work station 2, 
computational server 3, and database server 4. Two-way 

25 communication can occur between the servers such that the 
application server 1 receives input from the other 
servers and sends output information to the other 
servers. A user can interact with the system through a 
client workstation 2 which communicates with the 

30 application server, for example, by sending a query or 
command and by receiving the results of a computer 
implemented process of the invention. 



WO 03/106998 PCT/US03/18838 



27 

An application server 1 can extract data from 
the database server 4 or can launch simulations 
calculated on the computational server 3, for example, in 
response to a query or command received from the client 
5 workstation. Examples of- databases that can be accessed 
by the database server include a compound database, gene 
database, reaction database, bioparticle database or a 
reference database, each of which is described in further 
detail below. Simulations that can be accessed by a 
10 computational server 3 can include, for example, a single 
optimization analysis, deletion analysis, robustness 
analysis, phase plane analysis or time-course analysis 
each of which is set forth in further detail below. 

A multi-server architecture allows for the 
ability to manage information by storing the information 
on separate servers that can reside in the same location 
or can be globally distributed as in an application 
service provider (ASP) distribution model. The 
architecture can include any of a number of compatible 
network systems known in the art such as a local area 
network (LAN) or a wide area network (WAN) . 
Client-server environments, database servers and networks 
that can be used in the invention are well known in the 
art. For example, the database server can run on an 
operating system such as UNIX, where the operating system 
is running a relational database management system, a 
World Wide Web application or a World Wide Web server. 

Instructions or software code to implement a 
process of the invention can be written in any known 
30 computer language including, for example, an object 
oriented language such as Java or C++, a visual 
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programming language such as Visual Basic or Visual C++, 
or other languages such as C, FORTRAN or COBOL and 
compiled using any well-known compatible compiler. 

The software of the invention can be run from 
5 instructions stored or active in a memory, such as random 
access memory, on a host computer system. Similarly, 
information utilized in model construction and use, such 
as network components and network models/ is stored in a 
memory on a host computer system such as a read only 
10 memory. A memory or computer readable medium can be a 
hard disk, floppy disc, compact disc, magneto-optical 
disc, Random Access Memory, Read Only Memory or Flash 
Memory. A computer system that contains the memory or 
computer readable medium used in the invention can be a 
15 single computer or multiple computers distributed in a 
network. 

A database or data structure of the invention 
can be represented in a markup language format including, 
for example, Standard Generalized Markup Language (SGML), 
Hypertext markup language (HTML) or Extensible Markup 
language (XML) . Markup languages can be used to tag the 
information stored in a database or data structure of the 
invention, thereby providing convenient annotation and 
transfer of data between databases and data structures. 
In particular, an XML format can be useful for 
structuring the data representation of reactions, 
reactants and their annotations; for exchanging database 
contents, for example, over a network or internet; for 
updating individual elements using the document object 
model; or for providing differential access to multiple 
users for different information content of a data base or 
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data structure of the invention. XML programming methods 
and editors for writing XML code are known in the art as 
described, for example, in Ray, ^Learning XML" O'Reilly 
and Associates, Sebastopol, CA (2001) . 

5 The system architecture of Figure 1 is 

exemplary. Those skilled in the art will recognize that 
a process of the invention can be implemented on any of a 
variety of compatible architectures. For example, ■ the 
functions carried out by the servers can be consolidated 

10 into fewer servers or, alternatively, different functions 
or modules, such as those set forth below, can be tiered 
into a greater number of servers if desired. Although a 
single client desktop 2 is shown in Figure 1, it will be 
understood that the system can be readily modified to a 

15 multi-user distributed application to support 

collaborative network model construction or simulation, 
for example, by including multiple client desk tops that 
access an application server 1. 

A computer implemented process of the invention 
20 performs specified manipulations of data or information 
in response to a command or set of commands given by a 
user. A computer implemented process of the invention 
can be carried out by a computer system that provides an 
interface for a user to interact with the process by 
25 means of at least one use-case. A user is someone or 
something that interacts with a computer system from 
outside of the system. A use-case is a sequence of 
actions that a system performs, usually in response to a 
user command or input, that yields an observable output 
30 or result that is of value to a particular user. 

Accordingly, a computer system of the invention can 
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include any of the hardware components and compatible 
software set forth above such that the system contains 
executable instructions to carry out the computer 
implemented processes and use-oases set forth below. 

5 A use-case can be used to access or utilize a 

browser. A browser ' is understood to be a program which 
gives some means of viewing the contents of a data 
element in one or more database and of navigating from 
one data element to another. A data element can contain 
10 information about a compound, reaction, or organism and 
can be viewed, for example, by hypertext links accessed 
by the browser. 

An overview of an exemplary model construction 
process is provided in Figure 2 . Model construction is 

15 initiated 100 by selecting a bioparticle such as an 

organism, cell or virus or a biological system for which 
an in silico model is to be constructed. Although model 
construction will be described below with reference to a 
bioparticle for purposes of clarity, it will be 

20 understood that these steps can be carried out for a 
biological system within a bioparticle or encompasing 
more than one bioparticle. A bioparticle can be selected 
based on any of a variety of factors including, for 
example, the identification that it is a pathogen and the 

25 desire to create an in silico model for determination of 
effective therapeutic approaches to preventing 
pathogenecity, the identification that it is useful in an 
industrial process and the desire to create an in silico 
model for determination of optimal growth or production 

30 properties, or the identification that it is involved in 
a disease and the desire to create an in silico model for 
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identification of therapeutic targets for treatment of 
the disease. Any virus, prokaryote, bacteria, archaea or 
eukaryote for which sequence and or biochemical 
information is available can be modeled according to the 
5 invention. Specific examples of bioparticles that can be 
simulated by the models and methods of the invention 
include Arabidopsis thaliana , Bacillus subtilis , Bos 
taurus f caenorhabditis elegans, Chlamydomonas reihardtil , 
Danio rerio, Dictyostelium discoideum, Drosophila 

10 melanogaster , Escherichia coli , hepatitis C -virus, 
Haemophilus influenzae, Helicobacter pylori, Homo 
sapiens, Mus mus cuius , Mycoplasma pneumoniae, Oryza 
sativa, Plasmodium falciparum, Pnemocystis carinii, 
Rattus nor-vegicus , Saccha romyc e s cerevisiae, 

15 Schizosaccharomyces pombe, Takifugu rubripes, Xenopus 
laevis or Zea mays, and the like. 

The construction process can include a step 200 
of model requisition. At this step, preliminary 
evaluation can be made to determine whether to proceed 
20 with creating a new model, or to use an existing model, 
if present, that can be modified. At this step or any 
time prior to or during the process, individuals can be 
designated to have access to the model or the databases 
associated with the model can be selected. 

25 Access can be based on a particular set of 

rights provided to a user or set of users. For example, 
rights can include or exclude the ability to view all or 
part of the information stored in a database, the ability 
to edit all or part of the information stored in a 

30 database, . the ability to copy all or part of the 

information stored in a database, the ability to delete 
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all or part of the information stored in a database, the 
ability to use all or part of the use-cases included in a 
computer system, or a combination of these abilities. 
Limited access, for example, with respect to the right to 
5 edit stored information, can provide quality assurance 
and quality control of a database -and the information 
stored therein. Security and limited access rights can 
be achieved using known computer security algorithms and 
hardware such as those available from the SANS (System 

10 administration, networking and security) Insititute 
(available on the world wide web at sans.org) or 
Pentasafe (Houston TX, available on the world wide web at 
pentasafe.com). One or more users can be allowed access 
at a status of curator thereby having full rights 

15 necessary to access and maintain algorithms, models or 
databases . 

As shown in Figure 2, the model construction 
process can include a step 300 of collecting relevant 
organism specific information. At this step, a user such 

20 as a model developer can create a file structure for the 
bioparticle under which information relevant to the 
bioparticle is indexed and stored. Information that can 
be stored at this step includes, for example, a general 
description of the bioparticle, an appropriate taxonomy 

25 identification for the bioparticle that allows cross 
reference to information in databases or scientific 
publications or links to the NCBI Taxonomy Database 
(available on the world wide web at 
ncbi . nlm . nih . gov/ Taxonomy /taxonomyhome . html/ ) . 

30 At this step, a list of genes that encode gene 

products that perform reactions carried out by one or 
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more bioparticles of interest, for example, can be 
created. Many of these reactions occur due to the 
activity of a biomolecule catalyst or transporter, which 
are created through transcription and translation of the 
5 open reading frames (ORF) or genes found within the 
genome of a bioparticle. For purposes of brevity, 
reactions that occur due to the activity of a gene 
product and for which a cognate ORF is associated are 
referred to as gene-encoded reactions. Other reactions 

10 occur either spontaneously, through non-enzymatic 

processes or through proteins for which an ORF has not 
been associated are referred to as non gene-encoded 
reactions. Management of the data, for example, using a 
universal data management module can be achieved as 

15 described in further detail below. 

Every reaction whether or not it is gene- 
encoded contains one or many reactants, which are the 
chemical species or compounds involved in the reaction. 
These reactants can be designated as either substrates or 

20 products each with a discrete stoichiometrio coefficient 
assigned to them to describe the chemical con-versxon 
taking place in the reaction. The reactants are further 
specified according to the cellular compartments in which 
they are present. For example, in a reaction database, a 

25 distinction is made between glucose in the extracellular 
compartment versus glucose in the cytosol. Additionally, 
reactants in the reaction database can be specified as 
primary or secondary metabolites to assist in visual 
representations of large networks of metabolic reactions. 

30 Each reaction is also described by the 

direction in which it can proceed with the choices being 
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either reversible or irreversible. If a reaction is 
reversible then it is possible for not only the 
substrates to be converted into products, but also for 
the products to be converted into the substrates. 
5 Whereas an irreversible reaction is constrained to 

proceed only in the direction that converts substrates 
into products. 

At step 300 data elements specifying 
information regarding gene or genomic sequences, or 

10 attributes thereof, for a bioparticle can be obtained 
from an available source by, for example, downloading 
from a database into a gene index. The information 
included in this index can be downloaded from a public or 
private database or from an internal bioinf ormatics 

15 support service. Examples of databases from which gene 
• or genome information can be downloaded include the 
databases described above and in Baxevanis, supra, 2002. 
Sequences and annotations for a bioparticle genome or for 
genome fragments such as genes can be imported and stored 

2 0 in a gene database. The gene index includes structural 
information such as nucleotide sequences and genome 
annotation. Genome annotation includes identification of 
the location of ORFs and identification of homologies to 
other known genes. This information can be used to 

25 determine the function of the associated gene product (s), 
which can then be linked to the appropriate reactions 
that are catalyzed by the gene product (s). 



Although it is possible to access sequence data 
from outside databases during model construction and use, 
30 a gene index provides the advantage of direct access to 
data that may be dispersed in multiple non~associated 
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databases and the advantage of uniform storage or 
handling of information- for efficient cross-referencing 
and access. The system can include an algorithm and 
software code for importing genome sequences with or 
5 without supporting annotations into a gene index . 
Importation can be manually activated by a model 
developer or other user who identifies an updated genome 
dataset and has rights to edit a genome database or gene 
index. Alternatively, an algorithm and its implementing 

10 code can be included that automatically updates the 

information in a gene index by downloading information 
from an external database at a particular time interval 
or in response to a signal from the external database or 
its administrator that the data has been updated or 

15 modified. 



Also at step 300 other relevant information 
such as that available from the scientific literature 
regarding the genetics, biochemistry, cell biology and 
physiology of a bioparticle of interest can be gathered. 

20 These sources of information can be indexed in a citation 
library. The information is gathered in preparation for 
the process of constructing a network model which is 
described in detail below. The citation library can be 
integrated into a computer system that is used to make 

2 5 and use a network model such that the information in the 
citation library can be accessed from cross-references or 
hypertext links to network model components such as 
genes, biomolecules , reactions and compounds. 



Other network reaction components can also be 
30 stored in one or more data bases and accessed in a 
computer implemented process of the invention. For 
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example, a compound database can be used to store 
information relevant to biological compounds and 
reactants including substrates and products of reactions 
can be identified from the compound data base. A 
5 database accessed in a process of the invention can be 
specific to a particular organism strain, organism, 
species, family, phylum or kingdom. Alternatively a data 
base can be a universal database that contains genes, 
reactions, compounds or other information that is not 
10 exclusive for any subset of biological organisms. Thus, 
a universal reaction database or universal compound 
database is provided and can be accessed in a process of 
the invention. 

Referring again to Figure 2, the process can 

15 include a step 400 of constructing a data structure of 

network reaction components. A computer implemented step 
can be invoked to form a data structure associating 
network gene components with network reaction components. 
Such associations establish a data set specifying a 

20 network model of connectivity between network reaction 

components. For example, an ORF of a bioparticle can be 
selected and its gene sequence or other attributes 
identified. Such ORF data elements, either individually 
or together, specify data elements or data sets of a 

25 network gene component. The gene component can be 

associated directly, or used to identify its encoded gene 
product as a corresponding network reaction component. 
Obtained or identified network reaction components and 
their associated attributes, such as the reactants, 

30 enzymes or proteins that carry out the reaction, or mRNA 
encoding the enzyme or protein, similarly constitute data 
elements or data sets that can be incorporated into a 
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network model by association with gene components. All 
other associated relationships and attributes of 
identified gene and reaction components can similarly be 
incorporated into the network model by similar 
5 association. Such associations of gene and reaction 
components define the connectivity of gene product 
production and the connectivity and flow of reactions 
components of a network model of the invention. 

As described further below, the process of 
association can be repeated for inclusion of additional 
network components until a sufficient number of 
components have been identified to specify a functional 
group of interconnected or interrelated network members. 
Component attributes such as activity, substrates, 
products, reactants and stoichiometry serve to 
automatically associate, by natural biochemical 
relationships, the individual network components into an 
interconnected functional model. The natural 
relationships formed can be modified, for example, by a 
developer or user pf a network model of the invention. 
Therefore, the process of identifying, including and 
associating network components into a model of the 
invention serves to define the connectivity and flow of 
components and activity within the boundaries of the 
model itself. 

Association of data elements or sets of network 
gene components with corresponding data elements or sets 
of network reaction components can be performed by any 
computational method well known to those skilled in the 
30 art. For example, the individual data elements that make 
up the resultant data set can be associated using 
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relational tables. Alternatively, data elements can be 
associated using, for example, functions such as 
indexing, pointing, querying and the like. Similarly, 
combinations of these and other structures or functions 
5 can similarly be employed to associate network components 
included in a model of the invention. Further, the data 
elements can be partitioned within a database based on 
related characteristics or attributes or stored randomly. 
Alternatively, different databases can be used- to store 
10 categorized or uncategorized data elements. Therefore, 

associations of network components can be accomplished by 
any electronic linkage, physical archival form or 
combinations thereof . 

A data structure that is formed by the computer 

15 implemented process of the invention can be any physical 
or logical relationship among reaction components that 
supports flux balance analysis. Briefly, the data set 
consisting of associated data elements can be directly 
employed as a data structure of the invention. For 

20 example, the associated data set can be accessed by query 
and response from, for example, designated servers or 
specified server functions, and the associated data 
elements invoked as a single data structure during 
application of a network model of the invention. 

25 Alternatively, such associations can be further 

manipulated into secondary forms that can be accessed and 
utilized in the computer implemented methods of the 
invention. Such secondary forms can be created by, for 
example, further indexing, partitioning or the creation 

30 of subfiles and substructures of the data elements. For 
example, some or all of the associated data elements 
describing gene and reaction components can be 
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consolidated into a single data set. Where less than all 
of the data elements describing network components of a 
model of the invention are consolidated, it can be 
beneficial to maintain the associations and relationships 
5 to the original data elements and data sets to provide a 
continuous link to all characteristics and attributes of 
any particular network component represented by a data 
element. Maintaining such links provides an advantage of 
invoking computational processes on data elements 
10 relevant to network model performance while allowing 
manipulation of input, optimization and output of all 
data elements of any network component or any specified 
subset thereof. 

Specific examples of associations that can be 
15 constructed of network gene components and network 

reaction components by the computer implemented processes 
of the invention are described further below and in 
Example I. Figures 3 and 4 described therein set forth 
exemplary data elements specifying network components of 
20 a network model of the invention and their associations 
in both object model and database schema forms. Figure 3 
shows an associated object model specifying the 
participating classes of network component data elements 
and • associations in a network model of a bioparticle. 
25 Figure 4 shows ah associated database schema specifying 
the participating tables of network component data 
elements and associations in a biochemical network of a 
bioparticle . 

As shown in Figures 3 and 4, the network 
30 components can be organized into tables such as a table 

for reaction, reactant, molecule, protein, peptide, model 
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reaction, model version or gene. Within each table is a 
collection of records for attributes of the network 
component. For each record the fields are populated by 
the information added during network model construction 
5 as described below. 



A record can contain an attribute that is 
- represented in any appropriate format known in the art 
including, for example, a string, integer, float, 
character or boolean expression. String records are used 

10 for records that will have fields representing 
descriptions such as those for official name, 
abbreviation, direction, notes and discriminator in the 
table for the reaction network component. Boolean 
records are used to represent attributes for which one of 

15 two values is descriptive including, for example, whether 
a reaction is a transformation, translocation, unknown 
enzyme class, unknown transporter class or simulation 
reaction in the table for the reaction network component. 
Integer records can be used to denote numerical values 

20 such as the 5 T coordinate, 3 T coordinate, gene length and 
protein length occurring in the gene table. Examples of 
records that are represented as a float are molecular 
weight in the molecule table and coefficient including, 
for example, kinetic constants or binding constants in 

25 the reactant table. 

Exemplary associations between network 
components are indicated in Figure 3 and Figure 4. The 
associations can be utilized during various stages of 
model construction. For example, for the construction of 
30 a gene-protein association the tables that participate 
include the Peptide table, PeptideProteinAssociation 
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table, PepPepProteinAssociation table and Protein table 
as shown in Figure 4. The classes that participate in 
creation of a gene-protein association include Peptide, 
PeptideProteinAssociation and Protein. As another 
5 example of constructing an association using the tables 
and classes shown in Figures 3 and 4, a protein-reaction 
association is constructed using the Protein class, 
ProteinReactionAssociation class and ModelReaction class 
and using a Protein table, ProteinReactionAssociation 
10 table, ProtProtReactionAssociation table and 
ModelReaction table . 

Although the invention has been exemplified 
above with respect to a relational database, one of skill 
in the art will appreciate that the concepts presented 
15 herein may be applied outside of the relational database 
system of operation. In particular, the concepts are 
applicable in any database environment including for 
example an object-oriented database, hierarchical 
database or network database. 

20 A data set specifying network component 

associations can be transformed into a mathematical 
description of the network system being constructed. For 
example, in the specific case of modeling biochemical 
networks of a bioparticle, biochemical reactions of the 

25 network model can be transformed into a set of linear 

algebraic equations and inequalities. An inequality sets 
a constraint on a reaction that specifies an upper or 
lower boundary for the reaction. ' A boundary can specify 
a minimum or maximum flow of mass, electrons or energy 

30 through a reaction or can specify directionality of a 
reaction. A boundary can be a constant value such as 
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zero, infinity, or a numerical value such as an integer. 
Alternatively, a boundary can be a- variable boundary 
value . 

The set of' equations and inequalities 
5 constitutes a mathematical description of the referenced 
network model. A data structure of mathematical 
equations can be further represented as a stoichiometric 
matrix S r with S being an m x n matrix where m 
corresponds to the number of reactants or metabolites and 

10 xx corresponds to the number of reactions taking place in 
the network. Each column in the matrix corresponds to a 
particular reaction n, each row corresponds to a 
particular reactant m, and each S m element corresponds to 
the stoichiometric coefficient of the reactant m in the 

15 reaction denoted n. 

A stoichiometric matrix provides a convenient 
format for representing and analyzing a network model 
because it can be readily manipulated and used to compute 
network properties, for example, by using linear 

20 programming or general convex analysis. A network model 
data structure can take on a variety of formats well 
known to those skilled in the art so long as it is 
capable of relating components and reactions in the 
manner exemplified above for a stoichiometric matrix and 

2 5 in a manner that can be manipulated to determine an 

activity of one or more reactions using methods such as 
those exemplified below. Other examples of network model 
data structures that are useful in the invention include 
a connected graph, list of chemical reactions or a table 

30 of reaction equations. Such a table of chemical 
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reactions can further be annotated with kinetic 
information about the chemical reactions and 
transformations. Kinetic information can be accessed and 
used to apply differential equations to a network model 
5 of the invention or the reaction components therein to 
integrate over time. 

An exemplary process for implementing step 400 
is shown in Figure 5 . The process is initiated at step 
410 and proceeds to step 414 where a model is created or 

10 loaded. If an open edition of a desired network model is 
not present or accessible in the computer system, the 
process can proceed to step 418 in which an open edition 
is created and can then proceed to step 420. An open 
edition of a network model is one that is being generated 

15 or under construction. After sufficient improvement to 
the model content and preliminary testing the model can 
be saved as a versioned model to capture the current 
content of the model as a basis for future simulation 
studies. A versioned model is saved such that a copy of 

2 0 the versioned model is archived and the content of the 

archived model is secured or not substantively modified. 
If at step 414 an open edition is present and accessible, 
then the network model can be loaded into, for example, a 
computer processor or memory at step 416 and the process 

25 can proceed to step 420. It will be understood that a 
versioned model can also be loaded at step 416, for 
example, in order to create an updated or modified 
version of the model so long as at least one copy of the 
versioned model is archived and the model once opened at 

30 step 416 is stored as an open model until being saved as 
a new version. 
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At step 420 gene associated reaction components 
are added to a data structure of network reaction 
components. An exemplary process for implementing step 
420 is shown in Figure 6. As the data structure is being 
5 built, appropriate associations for each reaction to one 
or more related proteins and one or more related genes is 
assigned. These associations capture the relationships 
between the genes and proteins as well as between 
proteins and reactions. In some cases one gene codes for 

10 one protein which then catalyzes one reaction. However, 
often there are multiple genes which are required to 
create a protein and often there are multiple reactions 
that can be carried out by one protein or multiple 
proteins that can carry out the same reaction. These 

15 "associations can be captured by boolean logic operators 
such as X> AND" or "OR". These associations can also be 
captured in an association diagram as set forth below in 
the context of a model construction module. 

The representation of these associations in a 
20 network model of the invention provides the advantage of 
readily visualizing and determining the implications of 
adding or eliminating model content at the genetic, 
protein or reaction level in the context of making a 
network model or funning a simulation with a network 
2 5 model. In general, each of the genes in the gene index 

is evaluated for inclusion in or exclusion from a network 
model. If a gene is excluded, a reason can be provided 
in the annotations associated with the network model. 



The associations of network gene and reaction 
3 0 components can be implemented in a variety of different 

procedures. For example, the associations can be made in 
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or series. Additionally, a number of intermediate steps 
or groupings in the associations also) can be performed to 
facilitate or organize the resultant data structure. A 
5 specific example of the process of step 420 is where the 
implementing instructions invoke the selection or 
identification of a network reaction component based on 
an identified gene component. 

Identified gene components can be obtained, for 
10 example, from accessing a source of open reading frames 
(ORF) . The source can be derived from a variety of 
different resources and will depend on the network model 
intended to be constructed. For example, where a network 
model representing a biochemical pathway or a bioparticle 
15 function is to be constructed, a source of ORF data 
representing the activities of the pathway or the 
bioparticle function can be used. Specific examples 
include a gene database for the glycolysis pathway or a 
gene database for cellular metabolism. Similarly, where 
20 a network model representing the functions and activities 
of a bioparticle or subsystems thereof, a genomic 
database representing a substantially complete catalog of 
the bioparticle encoded genes can be used. 

One advantage of using an annotated network set 
25 of ORFs in constructing a network model of the invention 
is that it serves as an internal check on both the 
incorporation of network reaction components and on the 
completeness of the resultant model. For example, 
proceeding through a closed or finite list of gene 
30 components to be incorporated into a network model serves 
to internally constrain the number of possible 
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associations as well as identify aberrantly included or 
aberrantly associated network components. Therefore, 
construction of a network model from an annotated network 
set of ORFs provides both an upper and a lower limit for 
5 the components to be associated in the resulting data 
structure. Accordingly, model construction can proceed 
in a finite space of components and associations. 

Such a closed list of network gene components 
can be, for example, small such as would be for a pathway 

10 or bioparticle function. A closed list also can be, for 
example, large such as a bioparticle or organism genome. 
It is not necessary that an annotated network set be 
specified in a single list or file, or stored as a unique 
data entity. Instead, an annotated network set can be, 

15 for example, a subset of a larger database. Therefore, 
all that is required is the delineation of those ORFs 
included in an annotated network set from those excluded 
from the set . 

Regardless of the actual size of an annotated 
2 0 network set of ORFs, such a gene component set provides a 
genetic catalog or checklist for which the computer 
implemented process can proceed through and ensure that 
the listed gene components have been accounted for by, 
for example, either inclusion or exclusion from the 
25 network model being constructed. Additionally, the 

genetic catalog also can be used to invoke additional 
queries that call or proceed through routines relating to 
the identification and association of interactive and 
interrelated gene and reaction components. Invoking such 
30 routines or other analyses provides for a more complete 
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or thorough representation of the authentic system is 
reproduced in the constructed network model. 

For example, starting with a single ORF, the 
process of the invention can generate queries for 
5 identifying the corresponding encoded gene product and 
. attributes as well as any associated subunit components, 
their cognate ORFs and additional reaction constituents 
such .as substrates, products and cofactors." From that 
initial ORF and its identified gene products, cognate 

10 gene and gene product components, additional queries can 
be further invoked that expand on these associations by 
identifying network components related to the component 
being analyzed. Such expanded relationships can be, for 
example, the search and identification of network 

15 components upstream or downstream from the analyzed 

activity or physical interaction or of components and 
activities that are required to produce or deplete 
reaction constituents for the analyzed activity. 

Higher, "levels of expansion based on the 
20 initially selected ORF and its associated gene product 

can additionally be invoked depending on the need of the 
user or until queries and searches are exhausted. The 
computer implemented process can then proceed, for 
example, to the next ORF within the annotated network set 
25 to invoke the above queries and routines for 

identification of further reaction components and 
association into a network model data structure. 
Repeating this process of selecting an ORF, identifying 
its corresponding reaction component, querying and 
30 identifying interactive and interrelated gene, cognate 
gene and reaction components as well as reaction 
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is analyzed will yield a comprehensive group of network 
components that can be included by association into the 
network model being constructed. 

5 An additional advantage of model construction 

from an annotated network set of ORFs is that it provides 
or allows for the creation of data structure associating 
gene components with reaction components that will 
capture the inherent complexity of biochemical systems or 

10 living bioparticles . Moreover, such complexity can be ' 
reproduced in a network model with minimal knowledge or 
empirical determinations of the complete interactions or 
interrelationships of the system. Capturing and 
reproducing the natural complexity of biological and 

15 biochemical systems allows for a more accurate 

reproduction of the natural system in the resultant 
network model. 

For example, there can be instances where a 
particular biochemical .function is redundantly encoded in 

20 a bioparticle' s or organism's genome. Redundancy can 
therefore result in different gene products exhibiting 
similar function being represented in the repertoire of 
gene products. However, inclusion in a model of only a 
single gene product or activity can produce inaccurate or 

25 incomplete predictions because modification or 

perturbation of that single gene product or activity will 
not account for the substitutability of similar functions 
being present in the natural bioparticle or organism. A 
specific example augmenting the predictability of a 

30 network model by capturing the natural complexity of a 
biochemical system through inclusion of associations 
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between network gene and reaction components is described 
below in Example I. Thus, entirely different phenotypes 
can be observed depending on whether component redundancy 
is accounted for in a model. Other examples benefitting 
5 model reproduction and predictability of the authentic 
system by the inclusion of gene component associations 
include, for example, characterization of epistatic 
effects, evaluation of regulation at the gene, protein 
and reaction levels, comparative evaluation of the 
10 activity of isozymes or determination of the completeness 
with which the subunits of a multimeric protein are 
present in a network model. 

Referring again to Figure 6, and with reference 
to the initial process of selecting and associating gene 

15 and reaction components within a data structure, the 
computer implemented process of the invention gathers 
information related to the selected ORF in an initial 
screening or triage step. This initial step focuses on 
identifying and including network components specific for 

20 the model desired to be constructed. The information can 
be gathered by, for example, querying the user, a 
database or a server and obtaining replies that yield in 
the alternative a decision to either include or exclude 
the selected gene component in the data structure. For 

25 example, positive answers to whether the gene function is 
known, it is within the scope of the model being 
constructed or to non-ambiguous annotation or gene 
attribute information allow for inclusion of the selected 
gene component into the developing model. In this 

3 0 regard, a gene component can have a known function and 
clear annotation of attributes but be outside the scope 
of the model and be excluded such as when a metabolic 
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model is being constructed but the selected ORF encodes a 
nucleic acid binding protein or vice versa. 

Once a network gene component is determined to 
be included within a model being constructed, the process 
5 queries the user or a data source for identification of 
its encoded gene product. Alternatively, the process can 
electronically translate the gene component nucleic acid 
sequence data and include that information directly, or 
search a gene product data base to obtain the encoded 

10 .amino acid sequence as well as other attributes. As a 

maintenance procedure of the system, those gene products 
not represented in the corresponding database can be 
deposited in the system at this point or marked for later 
deposit during routine maintenance procedures. Following 

15 identification or generation of the corresponding gene 
product information, the resulting gene and reaction 
components are associated into a data structure. 
Generally, such association can be accomplished by 
employing relational databases and tables. However, and 

20 as described previously, essentially any means known to 
those skilled in the art can be used to form such 
associations. 

Once a network reaction component is associated 
with a gene component, the process. can further implement 

25 the selection of a new ORF from the annotated network set 
of ORFs and proceed with identification of its encoded 
gene product and related attributes. The initial 
selection queries for determining inclusion or exclusion 
is performed as described above. Further, the selection 

30 of subsequent ORFs and their encoded gene products can be 
performed, for example, sequentially, in parallel or in 
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series with the previous or subsequent ORF selections and 
processing. The newly identified network reaction 
components can again be subsequently incorporated into 
the network model by association with its corresponding 
5 gene component. Additionally, the functional and 

characteristic attributes of the reaction components also 
can be incorporated into the data structure of the 
network model being constructed. 

As described previously, once a network 

10 reaction component is associated with a gene component, 
the process can proceed further to extract or query data 
repositories or the user for related gene and reaction 
components as well as associate attributes of the 
identified network reaction components. Such related 

15 components include identifying and associating, for 
example, functional activities such as biochemical 
reactions, binding properties and other functional 
attributes; reaction constituents such as reactants, 
products and cof actors; constituent gene products such as 

20 subunits and regulators, as well as the various network 
gene and reaction components for such additionally 
identified network components. The implementation of 
these routines also is shown in Figure 6. Finally, for 
each identified reaction component, the process of the 

2 5 invention additionally queries whether the gene product 
catalyzes or participates in other reactions or 
processes. This step serves to expand the model 
construction process at each component to higher levels 
of component search, identification and association. 

30 Therefore, for each ORF included in the model 

construction as a gene component, the computer 
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implemented process of the invention proceeds through 
routine 420 one or more times until responses to the 
decision points are negative or exhausted or until the 
productivity of the output is outweighed by burden on 
5 computer or user resources. The repetition of routine 
420 begins- at the square box in Figure 6 denoting 
inclusion of a gene in the model. Upon termination of 
routine 420 for a particular included gene component, the 
process of the invention can continue through the 

10 annotated network set of ORFs by selecting another ORF 

and subjecting it to the preliminary decision points for 
inclusion into the developing model. Once included as a 
gene component, routine 420 is again implemented to 
identify and associate its encoded gene product as a 

15 reaction component, cognate gene components, gene product 
subunits, reaction constituents, additional gene products 
participating in the identified activity and the like. 
The complete routine 420 process can be, for example, 
repeated one or more times until the constituent ORFs of 

20 annotated network set, or a functional subset thereof, 
are processed and analyzed in similar fashion. 

Therefore, the invention provides a data 
structure that can be formed in a process of the 
invention by the steps of (a) selecting an ORF from the 

25 annotated network set encoding a gene product having a 
network reaction function; (b) forming a data structure 
including the selected gene product, the data structure 
associating network gene components and network reaction 
components including cognate ORFs, encoded gene products, 

30 network reactions and reaction constituents, and (c) 

repeating steps (a) and (b) selecting another ORF from 
the annotated network set until substantially all of the 
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network gene components of the annotated network set have 
been surveyed for encoding a gene product having a 
network reaction function to produce a data structure 
establishing a data set specifying a network model of 
5 connectivity and flow. , The process can further include 
the steps of (a) determining the occurrence of a 
constituent gene product for the selected encoded gene 
product; (b) determining the occurrence of an additional 
gene product participating in the network reaction; (c) 
10 determining the occurrence of an alternative network 
reaction exhibited by a surveyed gene product; and (d) 
incorporating identified constituent gene products, 
participating gene products or alternative network 
reaction into the data structure. 

15 A process of the invention can further include 

a step of elemental balancing at least one network 
reaction. Similarly, a process of the invention can 
include a step of charge balancing at least one network 
reaction. Such balancing takes into account conservation 

20 of mass, elements and charge as they occur in a 

biological system. Upon entry of a reaction by a user, a 
routine can be implemented to compare the substrates and 
products of a reaction to determine if mass is balanced 
such that the number of each atom type that enters a 

25 reaction in the substrates, matches the number that exits 
the reaction in the products. A similar comparison of 
the charge on substrates and products can be used to 
automatically determine if charge is balanced in a 
reaction that has been entered into the network model. 

30 If charge and mass are balanced the process is allowed to 
proceed to the next step in the construction process. 
However, if imbalance is found the system can send an 
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appropriate message to the user indicating that the 
reaction is not balanced. The message can further 
indicate the nature of the imbalance and suggest reaction 
constituents to add or remove in order to satisfy mass or 
5 charge balance. 

By monitoring the balance of charge, elements 
and mass on the reaction network the system makes 
resources available to a user that allow the user to 
interactively construct a network model that reflects the 
flux of mass and charge in a biochemical reaction network 
or biological system. Although mass, elements and charge 
balancing is not necessary for all applications of the 
network models of the invention, establishing this 
balance can account for phenotypes or system behaviors 
that occur in response to the net consumption or 
production of charge or a particular element. For 
example, the production of protons can affect cellular 
processes by altering pH, changing membrane potential, or 
contributing to processes that are energetically effected 
by proton influx/efflux such as metabolite transport and 
ATP levels. 

The process of constructing a data structure of 
network reaction components can include a step 430 of 
incorporating a network reaction that is not gene-encoded 
25 and corresponding reaction constituents into a data 
structure of network reaction components as shown in 
Figure 5. While many of the reactions of a bioparticle 
are associated with genes, there can also be a number of 
reactions included in a model for which there are no 
30 known genetic associations. A non gene-encoded reaction 
can be identified, for example, from the biochemical 
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literature or identified during the course of model 
construction based on the need for a reaction to satisfy 
a macro requirement deficiency. Knowledge of a gene or 
biomolecule that is associated with a reaction in a 
5 network model of the invention is not required fpr 

simulation using the model. However , such information 
provides advantages for efficient model building and for 
evaluating the results of a simulation. 

At step 430 reactions that occur spontaneously, 
10 that are not protein-enabled or that have not been 

associated with a particular gene product or open reading 
frame can be added to a data structure of network 
reaction components. Alternatively, a reaction can be 
added absent biological evidence indicating the 
15 occurrence of the reaction in a system being modeled, for 
example, based on results of a simulation and the 
identification of the need to satisfy a macro requirement 
deficiency by adding the reaction. 

One or more non gene-encoded reactions can be 
2 0 added to a network model during the course of model 

construction. Such a reaction can be associated with 
other reaction components such as reaction constituents 
and, where known, a cognate protein. The process can be 
carried out in the context of the model content browser. 
25 The computer implemented process is initiated when a 

determination is made to add a non gene-encoded reaction 
to a reaction index. The determination can be made by 
querying a user and obtaining a reply that yields an 
alternative decision that the reaction does or does not , 
30 exist in a reaction database. If the reaction occurs in 
a reaction database to which the user has been given 
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access, the reaction can be selected by the user and the 
system will automatically include the reaction in the 
reaction index. Alternatively, if the reaction does not 
exist in the reaction database, the user can be queried 
5 to enter the reaction and its corresponding reaction 
constituents into the reaction index. 

A reaction that is added to a reaction index 
can be added to a reaction database. The system can be 
configured to automatically add the reaction to the 
reaction database. Alternatively, the reaction can be 
displayed to a curator who responds to a query regarding 
whether or not the reaction is to be added to the 
reaction database. If the curator responds in the 
affirmative, the computer implemented process can add the 
reaction to the reaction database. Alternatively, a 
negative response by the curator will prevent addition of 
the reaction to the reaction database at that time. The 
process can proceed to query the user to edit reaction 
details such as the confidence level or to add a 
reference citation. 

The reactions in a data structure of network 
reaction components can be assigned to subsystems if 
desired. The use of subsystems provides advantages for a 
number of analysis methods such as pathway analysis and 
25 can make the management of model content more efficient. 
The model developer can specify the name of a subsystem 
and then assign reactions to the subsystem. This 
assignment allows a user to search for reactions in a 
particular subsystem which may be useful in performing 
30 various types of analyses. Furthermore, assignments of 
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subsystems can be indicated on reaction maps, thereby 
facilitating evaluation of simulation results. 

The reactions included in a data structure of 
network reaction components can be obtained from a 
5 reaction database using use-cases that are, for example, 
set forth below. Alternatively, reactions can be newly 
added, for example, by obtaining compounds from a 
compound database and building a reaction using methods 
similar to those set forth above for creating a reaction 
10 database. Reactions added at this stage of model 

construction can be subsequently added to a reaction 
database . 

The reactions added in steps 420 and 430 are 
intra-system reactions. Intra-system reactions are the 

15 chemically and electrically balanced interconversions of 
chemical species and biochemical processes, which serve 
to replenish or drain the relative amounts of certain 
metabolites. These intra-system reactions can be 
classified, for example, as either being transformations 

20 or translocations. A transformation is a reaction that 
contains distinct sets of compounds as substrates and 
products, while a translocation contains reactants 
located in different compartments. Thus, a reaction that 
transports a metabolite from the extracellular 

25 environment to the cytosol, without changing its chemical 
composition is classified as a translocation, while a 
reaction such as the phosphotransferase system (PTS) 
which takes extracellular glucose and converts it into 
cytosolic glucose- 6-phosphate is a translocation and a 

30 transformation . 
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Referring again to Figure 5, the process of 
constructing a data structure of network reaction 
components can include a step 440 of incorporating an 
exchange reaction for an external reaction component and 
5 corresponding reaction constituents into a data 

structure. Exchange reactions are the reactions that 
will allow compounds to be introduced and removed from 
the network for the purposes of simulation. Exchange 
reactions can be created based on empirically observed 
10 phenotype or behavior of a biological system. 

The metabolic or other biochemical demands 
placed on a biological system can be readily determined 
from the dry weight composition of a cell which is 
available in the published literature or which can be 

15 determined experimentally. The uptake rates and 
maintenance requirements for an organism can be 
determined by experiments in which the uptake rate is 
determined by measuring the depletion of the substrate 
from the growth medium. The measurement of the biomass 

20 at each point can also be determined, in order to 
determine the uptake rate per unit biomass. The 
maintenance requirements can be determined from a 
chemostat experiment. For example, the glucose uptake 
' rate can be plotted versus the growth rate, and the 

25 y-intercept interpreted as the non-growth associated 
" maintenance requirements. The growth associated 

maintenance requirements are determined by fitting the 
model results to the experimentally determined points in 
the growth rate versus glucose uptake rate plot. A data 

30 set of the invention can be modified to enumerate these 
experimentally determined demands using exchange 
reactions . 
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Exchange reactions are those which constitute 
sources and sinks, allowing the passage of metabolites or 
other network components into and out of a compartment or 
across a hypothetical system boundary. These reactions 
5 are included in a model for simulation purposes and 

represent the metabolic demands placed on an organism. 
While they may be chemically balanced in certain cases, 
they are typically not balanced and often have only a 
single substrate or product. As a matter of convention 
10 the exchange reactions are further classified into demand 
exchange and input/output exchange reactions. 

i 

Step 440 of a computer implemented process of 
the invention can be carried out in an exchange reaction 
browser. The computer implemented process can include a 

15 routine where input/output exchange reactions are added 
for extracellular reactants. The extracellular reactants 
in the data structure can be automatically displayed on a 
graphical user interface based on their identification 
during steps 420 and 430. The process can proceed to 

20 query the user whether or not to add input/output 
exchange reactions for all reactants that are 
extracellular. If the user answers in the affirmative, 
the process proceeds to insert exchange reactions for all 
extracellular reactants. Alternatively, if the user 

25 answers in the negative, the user is given access to 
evaluate the extracellular reactants and is further 
queried as to whether each should have an input/output 
reaction added. 



Thus, for each of the extracellular metabolites 
30 a user can specify or create a corresponding input or 
output exchange reaction. Generally, the system will 
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represent these reactions as reversible with the 
metabolite indicated as a substrate, a stoichiometric 
coefficient of one and no products produced by the 
reaction. This particular convention is adopted to allow 
5 the reaction to take on a positive flux value for its 
activity level when the metabolite is being produced or 
drained out of the system and a negative flux value when 
the metabolite is being consumed or introduced into the 
system. These reactions can be further constrained 
10 during the course of a simulation to specify which 

metabolites are available to the cell and which can be 
secreted by the cell. 

A demand exchange reaction can be introduced 
for any reactant in a network model of the invention. 
These reactions are introduced for biochemical demand 
constituents which are reactants that are required to be 
produced by the cell for the purposes of creating a new 
cell such as amino acids, nucleotides, phospholipids, and 
other biomass constituents, or metabolites that are to be 
produced for alternative purposes. A demand exchange 
reaction is generally specified as an irreversible 
reaction containing at least one substrate. These 
reactions are typically formulated to represent the 
production of an intracellular component by the metabolic 
network or the aggregate production of many reactants in 
balanced ratios such as in the representation of a 
reaction that leads to biomass formation, also referred 
to as growth. 

At step 440, the computer implemented process 
30 can also include a routine where demand exchange 

reactions are added for biomass constituents. The 
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process can proceed to query the user whether or not to 
add demand exchange reactions for all reactants that are 
biomass constituents. If the user answers in the 
affirmative, the process proceeds to insert demand 
5 exchange reactions for all biomass constituents. 

Alternatively, if the user answers in the negative, the 
user is given access to evaluate the biomass constituents 
and is further queried as to whether each should have a 
demand exchange reaction added. 

10 Generally, the system will represent these 

reactions as irreversible and specify the reactant as a 
substrate with a stoichiometric coefficient of unity. 
With these specifications, if the reaction is active it 
leads to the net production of the reactant by the 

15 network model due to potential production demands. 

Examples of processes that can be represented as a demand 
exchange reaction in a network model data structure and 
analyzed by the methods of the invention include, for 
example, production or secretion of an individual 

20 protein; production or secretion of an individual 

metabolite such as an amino acid, vitamin, nucleoside, 
antibiotic or surfactant; production of ATP for 
extraneous energy requiring processes such as locomotion; 
or formation of biomass constituents. 

25 The process of constructing a data structure of 

network reaction components can include a step 450 of 
creating one or more aggregate demand exchange reactions, 
which specify an aggregate reactant demand flux. 
Aggregate demand exchange reactions are demand exchange 

30 reactions that utilize multiple reactants in defined 

stoichiometric ratios. An example of an aggregate demand 
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reaction is a reaction used to simulate the concurrent 
growth demands of production requirements associated with 
cell growth that are placed on a cell, for example, by 
simulating the formation of multiple biomass constituents 
5 simultaneously at a particular cellular growth rate. 
Thus, an aggregate reactant demand flux can define a 
phenotypic output for growth. Other phenotypic outputs 
that can be defined by an aggregate reactant demand flux 
include, for example, biomass production, energy 

10 production, redox equivalent production, catabolite 

production, biomass precursors, polypeptide production, 
amino acid production, purine production, pyrimidine 
production, lipid production, fatty acid production, 
cofactor production, production of a cell wall component 

15 or transport of a metabolite. 

Step 450 of a computer implemented process, in 
which aggregate demand exchange reactions are 
constructed, can be carried out in an exchange reaction 
browser. A routine can be implemented in which the 

20 -reactants in the reaction database are automatically 
displayed on a graphical user interface. A user can 
review the contents of the display and identify reactants 
to be included in an aggregate demand exchange reaction. 
Biomass demand exchange reactions can be sequentially 

25 added to the aggregate reaction and biomass constituents 
can be added to the aggregate reaction. The user can be 
queried as to whether additional reactants should be 
added to the reaction. If the response is in the 
affirmative, additional reactants can be added. 

30 Alternatively, if the response is negative, the computer 
implemented process can specify stoichiometric 
coefficients for all reaction participants. The user can 
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then be queried to add additional aggregate exchange 
reactions. The user can repeat the process from the step 
of adding additional biomass demand exchange reactions. 
The routine can be repeated until a desired number of 
5 aggregate demand exchange reactions have been added. 

Therefore, the invention provides a computer 
implemented process for constructing a scalable output 
network model of a bioparticle. The process includes the 
computer implemented steps of: (a) accessing a database 

10 of network gene components including an annotated network 
set of open reading frames (ORFs)of a bioparticle genome; 
(b) forming a data structure associating the network gene 
components with network reaction components, the data 
structure establishing a data set specifying a network 

15 model of connectivity and flow of the network reaction 
components; (c) modifying the data set to enumerate a 
biochemical demand on the specified network model, and 
(d) transforming the modified data set into a 
mathematical description of reactant fluxes defining the 

20 network model of connectivity and flow, wherein the 

enumerated biochemical demand corresponds to an aggregate 
reactant demand flux defining a phenotypic output of the 
network model of a bioparticle. 

Once intra-system and exchange reactions have 
25 been added to a data ■ structure of network reaction 

components, the process can move to step 460 in which 
testing is performed to identify network gaps or other 
macro requirement deficiencies. This primarily includes 
testing to locate gaps in the network or "dead-ends" in 
30 which a reactant can be produced but not consumed or 

where a reactant can be consumed but not produced. The 
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determination of these gaps can be readily calculated 
through the appropriate queries of a reaction index and 
need not require the use of simulation strategies, 
however, simulation analyses are a possible approach to 
5 locating such metabolites. Gaps in a reaction network 
model can be identified by examining each of the 
reactants in the model to determine if they can be 
consumed and produced by the reactions therein. Gap 
analysis is accomplished using an algorithm that 

10 determines for each reactant if it occurs only once as a 
reactant or occurs multiple times as only a substrate or 
product when all the reactions are irreversible. If 
either of these criteria is satisfied then the reactant 
is displayed to a graphical user interface as a macro 

15 requirement deficiency. The user is then queried as to 
whether the gap should be accepted. The user can then 
decide to add or remove a reaction component from the 
network to eliminate the macro requirement deficiency, 
thereby incorporating an ameliorating network reaction 

2 0 component. Alternatively, the user can leave the macro 

requirement deficiency in the network if it is determined 
to have an insignificant effect on a simulation that is 
to be run using the network model or if the effects of 
the deficiency are to be determined in a simulation. 

25 An ameliorating network reaction component that 

is capable of augmenting competence of the connectivity 
and flow of a network model can be identified by a user 
that interacts with the network model in a computer 
implemented process, as set forth above. A computer 

30 implemented process can also identify the ameliorating 
network reaction component automatically. Thus, an 
algorithm that identifies a macro requirement deficiency 
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can further query a user to select , from a list of 
candidate reaction components, one or more reaction 
components that satisfy the deficiency. In the case 
where a macro requirement deficiency results in a 
5 reactant that is produced but not consumed, reactions 
from the universal reaction database that consume the 
reactant can be suggested as candidate ameliorating 
network reaction components. Alternatively, in the case 
where the macro requirement deficiency results in a 
10 reactant that is consumed but not produced, reactions 
from the universal reaction database that produce the 
reactant can be suggested as candidate ameliorating 
network reaction components. 

Alternatively, the computer implemented process 

15 can incorporate the ameliorating network reaction 

component automatically. Automatic incorporation can be 
achieved by an iterative process in which a candidate 
reaction component is tested in the network model, a gap 
analysis is performed and if the candidate reaction 

2 0 component augments competence of the connectivity and 
flow of the network model it is included or if the 
candidate reaction component does not augment competence 
of the connectivity and flow of the network model another 
candidate reaction is tested. The iterative process can 

25 be repeated until at least one reaction that augments 
competence of the connectivity and flow of the network 
model is identified. In the case that .more than one 
reaction is able to augment competence of the 
connectivity and flow of the network model, a user can be 

30 queried to make a selection or the selection can be made 
automatically based on criteria such as the confidence 
with which the reactions occur in other network models or 
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the presence of an ORF in a gene database that is 
-annotated to putatively encode one of the reactions. 

Thus, a process of the invention can include a 
5 step of incorporating an identified reaction component 
satisfying a macro requirement deficiency in structural 
architecture of a network model, wherein the 
incorporation supplements the connectivity and flow of 
the network model. For example, a process of the 

10 invention can include the steps of (a) determining the 
occurrence of a network reaction component satisfying a 
macro requirement deficiency in structural architecture 
of the network model, and (b) incorporating an identified 
network reaction component satisfying the macro 

15 requirement deficiency into the data structure to 

supplement the connectivity and flow of the network 
model . 

As shown in Figure 5, the process of 
constructing a data structure of network reaction 

20 components can include a step 470 of introducing 

confidence levels for reactions included in the data 
structure. The introduction of confidence levels 
enhances model specificity and provides the advantage of 
maintaining quality control and accountability for the 

25 content of the model. Accordingly, the reasons why a 
reaction is added or deleted from a model can be 
determined by the model developer contemporaneously, at a 
later date or by other users. Furthermore, a listing of 
evidence or reasons for including a reaction in a model 

30 can be maintained. 



WO 03/106998 



PCT/US03/18838 



67 

A step of annotating the reaction content of a 
model can be, for example, a dynamic activity that is 
ongoing throughout the model construction cycle and can 
be carried out at any stage of model construction. When 
5 a reaction is first added, a user such as the model 

developer can indicate the information levels and provide 
references. Alternatively, the user can add annotation 
details following entry of substantially all of the 
reactions to be included in a versioned model. 

In one embodiment, each reaction included in a 
data structure of network reaction components is 
annotated to reflect the confidence that the model 
developer has in the inclusion of the reaction in the 
model. The level of confidence is a function of the 
amount and form of supporting data that is available. 
This data can come in various forms including published 
literature, documented experimental results, or results 
of computational analyses. 



In the course of constructing a network model 
20 describing associations of network reaction components 

the types of data that will generally be accumulated and 
evaluated include, for example, biochemical data, genetic 
data, genomic data, physiological data, and modeling 
data. Biochemical .data includes information related to 
25 the experimental characterization of a chemical reaction, 
often directly indicating which biomolecule is associated 
with a reaction and the stoichiometry of the reaction or 
indirectly demonstrating the existence of a reaction 
occurring within a cellular extract. Genetic data 
30 includes information related to the experimental 

identification and genetic characterization of a gene 
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that encodes a particular biomolecule implicated in 
carrying out a biochemical event. Genomic data includes 
information related to the identification of an open 
reading frame and functional assignment, through 
5 computational sequence analysis, that is then linked to a 
biomolecule that performs a reaction. Physiological data 
includes information related to overall cellular 
physiology, fitness characteristics, substrate 
utilization, and phenotyping results, which provide 

10 evidence of the assimilation or dissimilation of a 
compound used to infer the presence of specific 
biochemical event including, for example, translocations. 
Modeling data includes information generated through the 
course of in silico modeling leading to predictions 

15 regarding the status of a reaction such as whether a 
reaction is needed to satisfy a macro requirement 
deficiency. 

The different forms of data elements that can 
be incorporated by association into a data structure of 

20 network reaction components, such as the data elements 
described above, can be ranked in terms of their 
importance toward determining the confidence level that 
will be assigned to a reaction. An exemplary ranking of 
highest information content to the lowest is as follows: 

25 biochemical, genetic, genomic, physiological, and 
modeling evidence . 

Within each type of data element or data set 
there are further hierarchies that can be established 
which can determine the overall quality of the data 
30 leading to an estimate that a particular form of data may 
provide no, low, medium, or high level of confidence. 
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Thus, confidence level can be determined from a 
hierarchical classification. Whether or not a reaction 
is included in a network model can be determined based on 
the relative confidence level in the hierarchy. For 
5 example, collectively hierarchical information levels can 
be used to heuristically determine an overall confidence 
level for a reaction in the model. A similar confidence 
scale could be used for other model content beyond just 
reactions. 

Depending upon whether or not information was 
gathered for each of the five relevant information types 
and, if information was gathered, the level of 
significance that the data holds with regard to the 
reaction, a score of no, low, medium, or high 
significance can be assigned. Additional annotation 
information in the form of textual notes can be attached 
to each reaction assignment as well as a list of relevant 
references gathered. Collectively these annotations, 
attached references, and the level of evidence associated 
with each of the data sources constitute the reaction 
rating details. 

A process of the invention can include a step 
of executing a heuristic logic decision algorithm that 
determines the level of confidence with which a network 
25 reaction component is included in a particular model. An 
overall reaction confidence level for the inclusion of a 
particular reaction in a data structure can be determined 
with a heuristic algorithm that evaluates the scores for 
information acquired in each of the five categories set 
3 0 forth above. In one embodiment, the overall confidence 
levels can range on a scale from one to five wherein 
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Level 1 means the reaction is speculative with no 
evidence,. Level 2 means the reaction is supported by 
minimal evidence , Level 3 means the reaction is supported 
by a fair amount of evidence, Level 4 means the reaction 
5 is highly probable with ample evidence and Level 5 means 
the reaction is certain to occur and has been validated. 
It is understood that these levels are exemplary and that 
a larger or smaller number of levels can be included to 
suit a particular application of the invention. An 
10 exemplary heuristic algorithm for determining confidence 
levels is described in Example II. 

These rating levels are provided as outputs 
such that they can be viewed by a model user or acted 
upon by a computational process when assessing the 

15' reaction content of a model. Thus, the confidence levels 
provide an annotation from which a model user can rapidly 
assess the confidence in a reaction assignment or 
identify groups of reactions listed at a particular 
confidence level. The user can be given access to 

20 investigate the reaction rating details if there is a 

need to further examine a particular reaction. In another 
embodiment, the level of confidence can provide a 
criteria for automatically determining inclusion or 
exclusion of a network reaction component in a network 

25 model. For example, a user can determine a threshold 
value such that reactions assigned ' greater confidence 
compared to the threshold value are automatically 
included in a network model while those reactions for 
which a lesser confidence level has been assigned are 

30 excluded from the model. 
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The process of constructing a data structure of 
network reaction components can include a step 480 in 
which a presimulation validation test is performed to 
determine if sufficient components of the network model 
5 are in place to allow simulation. A model validation 

report can be displayed to provide a general overview of 
the content of the model. The report can be reviewed 
before using the model for simulation and versioning. 
Examples of information that can be included in a 

10 validation report are ORFs that have been unevaluated for 
inclusion or exclusion from a model, ORFs included in the 
model that have "hypothetical", "unknown", or "none" 
included in their functional annotation, extracellular 
reactants that do not have an input/output exchange 

15 reaction included in the model or macro requirement 
deficiencies in the reaction network. Based on the 
displayed report a user can determine whether or not to 
modify an associated network model. 

A computer implemented process of the invention 
20 can further include a step of calculating a phenotypic 
output of a network model from its mathematical 
description. The phenotypic output can be calculated 
from the mathematical description using methods known in 
the art for flux balance analysis as described, for 
25 example, in Schilling et al . , J. Theor. Biol. 203:229-248 
(2000); Schilling et al . , Biotech. Bioeng. 71:286-306 
(2000); Schilling et al., Biotech. Proa. 15:288-295 
(1999) , and Varma and Palsson, Biotech. Bioeng. 
12:994-998 (1994). Briefly, a mathematical description 
30 such as a matrix or system of linear equations can be 

solved to calculate the null space that defines the set 
of steady-state metabolic flux distributions that do not 



WO 03/106998 



PCT/US03/18838 



72 

violate the mass, energy, or redox balance constraints. 
A point in this space represents a flux distribution and 
hence a phenotypic output for the network. An optimal 
solution within the set of all solutions can be 
5 determined using mathematical optimization methods when 
provided with a stated objective and a constraint set. 
The calculation of any solution constitutes a simulation 
of the model. 

The invention provides a computer implemented 

10 process for self-optimizing a network model of a 
bioparticle. The process includes the computer 
implemented steps of: (a) accessing a database of network 
gene components including an annotated network set of 
open reading frames (ORFs)of a bioparticle genome; (b) 

15 forming a data structure- associating the network gene 
components with network reaction components, the data 
structure establishing a data set specifying a network 
model of connectivity and flow of the network reaction 
components; (c) transforming the data set into a 

20 mathematical description of reactant fluxes defining the 
network model of connectivity and flow; (d) determining 
the competence of the connectivity and flow within the 
network model, the competence indicating underinclusion 
or overinclusion of network reaction component content of 

25 the network model, and (e) identifying an ameliorating 
network reaction component capable of augmenting the 
competence of the network model, incorporation of the 
ameliorating network reaction component into the data 
structure producing a modified data structure specifying 

30 in an optimized network model of the bioparticle. 
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Referring to Figure 2, the model construction 
process can include a step 500 of preliminary simulation 
testing and model content refinement. In this step the 
existing model can be subjected to a series of functional 
5 tests to determine if it can perform basic requirements 
such as the ability to produce the required biomass 
constituents and generate predictions concerning the 
basic physiological characteristics of the particular 
organism strain being modeled. Typically , the majority 
of the simulations used in this stage of construction 
will be single optimizations, which are set forth in 
greater detail below. Before a network model is used to 
examine the ability to use an aggregate demand reaction 
as an objective function, the model is typically tested 
to determine that it is capable of generating each of the 
individual components. As an example, before an 
aggregate flux to simulate growth is used, the model is 
examined to determine if all of the amino acids can be 
generated through the model reactions and inputs. Thus, 
the preliminary simulation testing involves the 
examination of the network to produce individual 
reactants by selecting the appropriate single demand 
exchange reactions as the objective and optimizing for 
the production of the reactant under a wide range of 
possible conditions. If the metabolite cannot be made 
then changes can be made to the model until a desired 
phenotypic characteristics such as growth can be 
simulated. 

Following a review of the content of the model 
and the results of preliminary simulation testing at step 
600 a decision can be made as to whether or not to 
version the network model. If the model is not 
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sufficiently complete to be versioned the process is 
repeated by returning .to step 500 or, if necessary 
another step in the process. Accordingly, model 
construction can be carried out in an iterative fashion 
5 in which steps of the process are repeated until a 

desired model is obtained. Once the network model is 
determined to be sufficiently complete the process 
proceeds to step 700 where the model is versioned. 
Iterative construction leads to the continuous 
10 improvement and refinement of in silico models. 

To make modifications to a model version a new 
open edition of the model can be created based on the 
model version that is to be modified. Once a model is 
versioned, it is generally not edited without creating a 
15 new edition. This includes changes to the reactions in 
the data structure of network reaction components and 
their associations to biomolecules and genes as well as 
changes to the reaction properties details such as the 
confidence level and references. 

2 0 The invention provides a system for 

constructing a scalable phenotypic output network model 
of a bioparticle. The system includes (a) an input data 
set of network gene components including an annotated 
network set of open reading frames (ORFs) of a 

25 bioparticle genome; (b) executable instructions forming a 
data structure associating the network gene components 
with network reaction components, the data structure 
establishing a data set specifying a network model of 
connectivity and flow of the network reaction components; 

30 (c) executable instructions modifying the data set to 

enumerate a biochemical demand on the specified network 
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model, and (d) executable instructions mathematically 
describing from the modified data set reactant fluxes 
defining the network model of connectivity and flow, 
wherein the enumerated biochemical demand corresponds to 
5 an aggregate reactant demand flux defining a phenotypic 
output of the network model of the bioparticle. 

The invention further provides a system for 
constructing a scalable phenotypic output network model 
of a bioparticle. The system includes (a) an input 

10 data set of network gene components including an 

annotated network set of open reading frames (ORFs) of a 
bioparticle genome; (b) executable instructions forming 
data structure associating the network gene components 
with network reaction components, the data structure 

15 establishing a data set specifying a network model of 

connectivity and flow of the network reaction components 
(c) executable instructions determining the occurrence o 
a reaction component satisfying a macro requirement 
deficiency in structural architecture of the network 

20 model, inclusion of an identified reaction component •■ 

satisfying the macro requirement deficiency in the data 
structure supplementing the connectivity and flow of the 
network model; (d) a heuristic logic decision algorithm 
determining confidence of the network reaction component 

25 within the data structure, and (e) executable 

instructions mathematically describing from the data set 
reactant fluxes defining the network model of 
connectivity and flow, wherein the mathematical 
description defines a scalable output network model of a 

30 bioparticle. 
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The invention provides a system for 
constructing a self-optimizing network model of a 
bioparticle. The system includes (a) an input data set 
of network gene components including an annotated network 
5 set of open reading frames (ORFs) of a bioparticle 
genome; (b) executable instructions forming a data 
structure associating the network gene components with 
network reaction components , the data structure 
establishing a data set specifying a network model, of 

10 connectivity and flow of the network reaction components; 
(c) executable instructions mathematically describing 
from the data set reactant fluxes defining the network 
model of connectivity and flow; (d) executable 
instructions computing competence of the connectivity and 

15 flow within the network model, the competence indicating 
underinclusion or overinclusion of network reaction 
component content of the network model, and (e) 
executable instructions augmenting the competence of the 
connectivity and flow within the network model, the 

20 executable instructions specifying inclusion or exclusion 
of an ameliorating network reaction component, wherein 
incorporation of the ameliorating network reaction 
component into the data structure produces a modified 
data structure specifying an optimized network model of 

25 the bioparticle. 

A computer system of the invention can include 
a number of separate modules that contain one or more 
use-cases having various functions associated with making 
and using 'a network model. One or more modules that can 
30 be included in the system include, for example, a 

universal data management module, model construction 
module, atlas management module, simulation module, data 
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mining, experimental data module, gene sequence analysis, 
module or any combination of these modules". A number of 
computer implemented processes of the invention are 
described below with reference to these modules. Those 
5 skilled in the art will understand that, although the 

modules provide particular advantages for organizing and 
managing information, as set forth below, the steps of a 
computer implemented process of the invention can be 
carried out with or without any or all of the modules. 

10 Network gene components can be stored in a gene 

index and partitioned into data elements and data sets 
each containing information identifying a particular gene 
with a name or genomic location and other information 
including, for example, structural information such as 

15 the primary ^sequence of the gene or annotations 

describing the structure or function of the gene. The 
data elements can be stored in such a way that when a 
network gene component is accessed or included in a data 
structure, information relevant to the gene is 

20 associated, for example, using a hyperlink. Thus, a step 
of accessing a database of network gene components can 
include accessing a network gene component and associated 
information stored in a particular data element. 

Information from which a network model is 
25 constructed or which can be used to modify an existing 
network model including, for example, a gene database, 
reaction database or compound database can be managed 
using a universal data management module. A universal 
data management module can include, for example, a use- 
30 case to maintain a citation library a use-case to 

maintain compounds, a use-case to maintain reactions, a 
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use-case to maintain bioparticle-specif ic data, or a 
combination of two or more of these use-cases . 

A use-case to maintain a citation library 
allows a user to manage references such as books, 
5 articles, journals and papers. This use-case can be 
performed using a third-party tool. The user can 
associate a reference with any particular reaction added 
to a model. - This use-case interacts with a user by 
providing the ability to add, delete, or edit any form of 

10 reference or citation that the user my wish t to include as 
part of a model for supporting information. The user 
enters a citation into the system, allowing the citation 
to be available for selection at any point when the user 
wishes to annotate any of the model content with a 

15 reference. 

A use-case to maintain a database such as a 
compound database, reaction database or bioparticle- 
specif ic database allows a user to access and edit data 
elements stored therein by adding, deleting or editing 

20 information relevant to a particular entry. Such use- 
cases interact with a user by displaying the contents of 
a database and allowing the user to add a new entry to 
the database, delete an entry from the database, or 
modify an existing entry. A modification of a compound 

25 database can include, for example, changing the atomic 

composition of a compound or adding, deleting or editing 
information such as physical properties listed in an 
entry for a particular compound. A modification of a 
reaction database can include, for example, changing the 

30 atomic composition of substrates and products, the type 
of reaction, stoichiometric coefficient for the reaction 
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or other information relevant to the reaction. A 
modification of a bioparticle-specif ic database can 
include, for example, changing names, taxonomic 
information, description of characteristic features or 
5 information regarding areas of practical application. A 
use-case for maintaining a database also provides a means 
to select a compound or reaction from a database, for 
example, using a command, query or index function that 
associates a selected compound or reaction to a network 
10 model data structure. 

A model construction module can be included in 
a computer system of the invention. The methods of the 
invention for constructing or generating a network model 
can be performed in a model construction module. This 

15 module provides use-cases for managing information 

regarding reaction content, properties of a biomolecule 
or set of biomolecules that catalyze a reaction, and 
nucleic acids encoding the biomolecules. The model 
construction module can be used for any stage of model 

20 construction and modification from initial assembly, to 
iterative model building, preliminary testing and 
versioning. A model construction module can .include, for 
example, a use-case to download a gene index, a use-case 
to maintain a gene index, a use-case to maintain model 

25 content, a use-case to evaluate a gene index, a use-case 
to maintain a reaction index, a use-case for model 
reconstruction, a use-case to maintain exchange 
reactions, a use-case to validate model structure and 
content, a use-case to gather model test data, a use-case 

30 to perform model testing, a use-case to version a model, 
a use-case to assign reactions to a region, or a 
combination of two or more of these use-cases. 



WO 03/106998 



PCT/US03/18838 



80 

A use-case to maintain model content allows a 
user to access and modify the content of model editions 
for a particular bioparticle or organism strain. This 
use-case interacts with a user by providing simultaneous 
5 access to a network model data structure, databases of 
relevant information and an association diagram. An 
association diagram is a display of associations between 
genes, the biomolecules they encode and reactions that 
are catalyzed or carried out by the biomolecules within a 
10 network model data structure. Exemplary association . 
diagrams are shown in Figure 7 . 

An association' diagram is updated in response 
to commands sent by a user to add, remove or otherwise 
modify the content of a network model data structure. 

15 Thus, the association diagram provides a convenient 

visualization tool for evaluating the effect of making 
changes at the gene, biomolecule or reaction level in a 
network model data structure. Take for example, a 
biomolecule catalyst having multiple subunits, where all 

2 0 of the subunits are required for activity "and. where each 
subunit is expressed from a different gene. Visual 
evaluation of the gene -biomolecule- react ion associations 
during model construction can allow a user to readily 
identify the full complement of genes required to perform 

25 a particular reaction. Thus, once any one of the genes 
is selected from the gene index for inclusion in a data 
structure the user can rapidly identify the full set of 
genes required to perform the reaction. Furthermore, 
because simultaneous access is provided to multiple 

30 databases, the identified information can be displayed to 
a user and the user can modify a data structure based on 
evaluation of the displayed information. 
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A use-case to maintain model content can also 
include commands to access and edit properties of a model 
edition such as its name, description and notes. The 
content of the model edition which can be viewed and 
5 modified includes the gene index, protein index, reaction 
index and associated references, exchange reactions, and 
network gaps. This use-case also provides algorithms to 
create a new model edition and change the properties of 
the edition- such .as .its name., ..description and notes. 

A gene index can be managed using a model 
construction module. A use-case to download a gene index 
allows a user to load into a computer system of the 
invention a gene index that has been generated from 
external third party software or downloaded from an 
external database. A gene index can be downloaded as a 
text file or in a spreadsheet and converted to a desired 
format using a suitable script. 

A use-case to maintain a gene index allows a 
user to access the data stored in a gene index 'and to 
edit the content of the data. This use-case interacts 
with a user by displaying the contents of a gene index 
and providing a means to, for example, modify the 
annotation and functional assignments made to individual 
open reading frames or genes within a genome. A gene can 
be added to a gene index or deleted from a gene index 
using this use-case. 

A use-case to evaluate a gene index allows a 
user to evaluate the gene index for a particular organism 
30 strain to determine the genes to be included in a model 
edition. This use-case interacts with the user by 
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displaying the contents of a gene index such that each 
gene or ORF can be evaluated for inclusion in a model 
edition. The user can send commands to eliminate a gene 
or ORF from the model or include a gene or ORF in the 
model. This use-case further prompts the user to 
indicate associations between genes , biomolecules and 
reactions . 

A use-case to maintain a reaction index allows 
a user to manage the reactions that are included in a 
model edition. This use-case interacts with the user by 
displaying the contents of a reaction index and providing 
a means to add a reaction to the reaction index, delete a 
reaction from the reaction index; add, remove or view a 
reference from a citation library associated with a 
reaction; assign a reaction to a subsystem; add a 
confidence level to a reaction, or annotate an entry for 
a reaction. 

A use-case for model reconstruction allows a 
user to determine the locations in a network model where 
0 a macro requirement deficiency or gap in the pathway 

structure occurs. This use-case interacts with the user 
by providing the ability to launch the gap analysis 
algorithm to locate reactants that are only consumed or 
produced in the network. The system then displays to the 
5 user a list of such metabolites along with information on 
whether they are only consumed or produced. The user can 
review and evaluate these macro requirement deficiencies 
and decide whether or not to take any action to eliminate 
the gap by addition or removal of reactions from the 
0 network. The user can iteratively add or delete 
reactions and rerun the gap analysis algorithm to 
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determine if the gap still exists. In addition the use- 
case can display candidate reactions that are potentially 
capable of satisfying an identified macro requirement 
deficiency. An exemplary process for identifying a macro 
5 requirement deficiency and adding a reaction component to 
satisfy the deficiency is provided in Example III. 

A use-case to maintain exchange reactions 
allows a user to manage the exchange reactions associated 
with a model edition. This use-case interacts with the 
10 user by providing access to a reaction index and allowing 
the user to identify reactions as an input exchange 
reaction, output exchange reaction or demand exchange 
reaction. In addition, a user can create, delete or 
modify an aggregate demand reaction with this use-case. 

Intra-system reactions can be managed with a 
use-case for maintaining model content while exchange 
reactions are managed by a separate use-case. 
Intra-system reaction components represent true 
biochemical reactions that occur in a bioparticle and are 
potentially associated with the genes in the bioparticle. 
Therefore, these reactions are subject to the assignment 
of associations between genes, proteins, and reactions. 
These reactions are typically atomically and electrically 
balanced. Additionally, confidence levels are only 
assigned for these reactions and not for exchange 
reactions . 

An algorithm can be included in a use-case for 
maintaining exchange reaction browser that automatically 
locates extracellular metabolites that occur in the N 
30 reactions that are included in a network model. 
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Extracellular metabolites identified by such an algorithm 
or any other means can be used for the creation of input 
or output exchange reactions. In addition, a use-case 
for maintaining exchange reactions can include an 
5 algorithm to locate biomass compounds or other 

biochemical demands and present them for the possible 
inclusion of biomass demand exchange reactions- The 
exchange reactions can be displayed such that a user can 
evaluate and select reactions to be included in a network 

10 model. Thus, the exchange reaction browser provides a 
means for a user to provide commands to exclude a 
reaction from a network model or to manually include a 
reaction that is not already present' in the universal 
reaction database. A reaction added to the network model 

15 will automatically be added to the reaction database and 
the reactants will be added to the compound database. 

A use-case to validate model structure and 
content allows a user to determine whether the structure 
and content of a model edition meet certain desired 

20 specifications before being versioned. This involves the 
completion of a number of basic structural analyses and 
the performance of some basic simulations to qualify a 
model as being valid. This use-case interacts with the 
user by performing a series of validation tests or 

2 5 queries on the contents of the model and reporting the 
results back to the user. The user can then view these 
results and if there are no significant problems 
identified, the model can be used for simulations and be 
versioned if desired. 

30 A use-case to perform model testing allows a 

user to refine the content of a model. In this stage the 
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existing model is subjected to a series of functional 
testing to determine if it can perform basic requirements 
such as the ability to produce the required biomass 
constituents and generate predictions concerning the 
5 basic physiological characteristics of the particular 

organism strain being modeled. A user interacts with this 
use-case by running simulations on the model. Based on 
the results of these simulations the user can make 
changes to the content of the model. 
10 Generally, the simulations used in this stage of 
construction are single optimizations. 

A use-case to version a model allows a user to 
version an open edition of a model. This use-case 
interacts with a user by saving an open edition of a 

15 network model as a versioned edition in response to 

commands given by the user. A versioned edition of a 
network model is saved such that no further changes can 
be made to the model version. A user assigned version 
number is given to each of the versions of a strain 

20 specific model. 

A use-case to assign or associate reactions 
relative to other components within a network model 
allows a user to identify a reaction as participating in 
a particular subset of reactions in a network such as in 

25 a particular metabolic pathway. The reactions in a 

network structure or reaction database can be subdivided, 
for example, according to biochemical or biological 
criteria, such as according to traditionally identified 
metabolic pathways (glycolysis, amino acid metabolism and 

30 the like) or according to mathematical or computational 
criteria that facilitate manipulation of a model that 
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incorporates or manipulates the reactions. Methods and 
criteria for subdividing a reaction database are 
described in further detail in Schilling et al . , J. 
Theor. Biol. 203:249-283 (2000). The use of subsystems 
5 can be advantageous for a number of analysis methods, 
such as extreme pathway analysis, and can make the 
management of model content easier. This use-case 
interacts with a .user by displaying the contents of a 
network model data structure and allowing the user to 
10 select a reaction and assign the selected reaction to a 
subsystem. 

A use-case to maintain constraint templates 
allows a user to maintain representative sets of data 
elements which define particular common intraparticle or 

15 environmental conditions. An example is a constraint 
template to represent aerobic growth conditions on 
glucose. A user interacts with this use-case by 
selecting a constraint template to be used as the 
baseline set of constraints used to run a simulation. 

20 The constraint template may be derived from a previous 
simulation as "well. This saves the user the time 
reguired to re-enter all of the constraints placed in a 
new simulation that was used for the same model in a 
previous simulation . 

25 Network model content also can be viewed or 

represented with maps that indicate the connectivity of 
reactions or fluxes that are present in the network. The 
maps can be output in a variety of different formats 
including, for example, two-, three- or multi-dimensional 

30 maps, diagrams and atlases/ Thus, the invention provides 
an algorithm for displaying a map of the reactions 
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included in all or part of a network model. A user can 
design a map by selecting reactions to be displayed on a 
map. Reactions are typically displayed with each of the 
reactants shown as nodes and the reactions connecting 
5' these reactants shown as arrows. The user can then 

arrange these reactions in a familiar layout on the map 
or can select to have the map layout automatically 
generated based on well established algorithms for 
auto-layout of graphs. Alternatively, an inverse map can 

10 also be designed wherein each of the reactions is 

indicated by a node while the metabolites are represented, 
by arrows connecting the two nodes. An inverse map is a 
different way to view a metabolic reaction network that . 
can offer advantages for the visualization of network 

15 function. 

A map can be further enhanced to show the flux 
of network components, biochemical demands, or aggregate 
demand through the reactions of a network based on the 
results of one or more simulation. Direction of flux can 

20 be represented by arrows or apparent directional movement 
of an image between reactants. The amount of flux 
through reactions of a network can be represented in a 
map, for example, by the relative width of reaction 
arrows where a gradient of arrow widths is correlated 

25 with the amount of flux, a color gradient correlating 

colors in a spectrum with the relative amount of flux or 
the rate at which apparent directional movements of an 
image occur between reactants. 



Also provided is a means for displaying a map 
30 that associates reactions with the biomolecules that 
carry out the reactions or the genes that encode the 
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bioraolecules . A map can further associate reactions, 
biomolecules and genes. 

An atlas management module can be included in a 
computer system of the invention and used to manage 
5 network maps and to organize them into a collection 
referred to as an atlas. An atlas is a collection of 
maps that can cover reactions spanning one or more 
organism. An atlas management module can contain a use- 
case to manage atlases and maps, a use-case to design a 
10 map, and a use-case to view and test a map. 

A use-case to manage atlases and maps allows a 
user to organize maps into atlases and allows the user to 
create or delete maps and atlases. This use-case 
interacts with a user by displaying a list of maps such 
15 that the user can add, delete or modify the collection of 
maps that are in a particular atlas. In addition, a user 
can interact with this use-case by copying an atlas, or 
map for efficient generation of a new map. 

A use-case to manage atlases and maps provides 
20 access to an atlas of maps contained in separate elements 
or folders within an atlas. Each bioparticle or organism 
strain can be correlated with a default map or set of 
maps so that when simulations are performed in a 
particular model, an appropriate map is first displayed. 
25 However , maps themselves need not be linked to models. 

Accordingly, a computer system of the invention provides 
a means to load any map and view any simulation result on 
the map, regardless of the organism (s) from which the map 
was generated. This functionality allows comparison of 
30 multiple simulation results from the same or different 
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models on the same map. Color scales can be used to 
represent different parameter values obtained from 
different simulations when displayed on the same map. 

A use-case to design a map allows a user to 
5 design maps of network models. These maps provide a 
convenient visual tool for evaluating the content of a 
model in terms of the reactions included in the model and 
how they are connected to one another. This, is a drawing 
and design tool that allows a user to design maps that 
10 represent network models at any of a variety of levels of 
detail from maps of individual pathways such as purine 
biosynthesis, to larger regions such as amino acid 
metabolism, and even substantially complete system maps 
of cellular metabolism. 

15 The design use-case interacts with a user by 

displaying a list of reactions included in a network 
model data structure and providing a canvas for graphic 
manipulation of map content. In response to a command 
from a user to include a reaction in a map, the use-case 

20 will automatically add the reaction to an appropriate 
location according to the connectivity of the network 
model data structure. The user can manipulate the map by 
altering the location of substrates and products and 
arrows connecting them will be redrawn consistent with 

25 the new location on the map and the connectivity of the 
network model data structure. Common data elements 
representing the same metabolite can be merged such that 
locations in the map where a particular metabolite occurs 
are connected or otherwise correlated or common elements 

30 can be kept separate on the map. Additionally, this use- 
case allows a user to send a command to render one or 
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more reactions that are present in a map as either 
visible or invisible. 

The design use-case can provide a user with 
analysis capabilities .to compare reactions placed on a 
5 map with reactions that occur in a particular model or 
region within a model. Visual features of the maps can 
include connectivity lines, options to handle secondary 
metabolites , hyperlinks to other maps, placeholders for 
numerical simulation results, or annotations. . Additional 

10 analysis features can be included on a map such as the 
ability to select a metabolite of interest and 
simultaneously view all of the reactions in which the 
metabolite participates. Analysis tools such as the 
visual features of the maps assist the user in 

15 determining the reactions which need to be placed in the 
map by providing a view of the connectivity of reactions 
in the network while allowing access to information 
databases that are useful in evaluating the properties of 
a particular reaction in the network. 

20 The maps can be used to display results from 

simulations and empirical data allowing for comparisons 
between simulations and experiments. For example, 
empirically determined results of gene expression, 
protein expression, protein-protein interactions or 

25 reaction rates can be compared to an in silico predicted 
flux distribution. 

Simulations can be performed and managed with a 
simulation module. This module contains use-cases for 
different types of simulations including, for example, 
30 single optimization, deletion analysis, robustness 



WO 03/106998 



PCT/US03/18838 



91 

analysis, phase plane analysis or time-course analysis. 
A simulation module can include, for example, a use-case 
to load or create a project, a use-case to manage 
simulations, a use-case to define optimization 
5 constraints, a use-case to perform a single optimization, 
a use-case to view single optimization results, a use- 
case to perform a deletion analysis, a use-case to view 
deletion analysis results, a use-case to perform 
robustness analysis, a use-case to view the results of 

10 robustness analysis, a use-case to perform phase plane 
analysis, a use-case to view results of phase plane 
analysis, a use-case to perform time-course analysis, a 
use-case to view results of time-course analysis, a use- 
case to compare simulation results, a use-case to compare 

15 single optimization and experimental results, a use-case 
to export simulation results or a combination of two or 
more of these use-cases . 

Simulations can be managed using use-cases to 
load/create, manage and export simulations respectively. 

20 A use-case to load/create a project allows a user to 

create scientific projects and assign them to a program. 
Each project can contain simulation studies and 
additional information that are related to a particular 
bioparticle or related to many bioparticles . Simulation 

25 studies contain the details of individual simulations and 
experiments. A use-case to load/create projects 
interacts with the user by displaying a list of available 
projects from which one or more can be selected and 
opened by the user. A user can organize and annotate 

30 simulation results or experimental data using a use-case 
to manage simulations. This use-case interacts with the 
user by allowing the user to edit the name of a project, 
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alter the program to which it belongs or annotate the 
project or program. A use-case to export simulation 
results can be used to convert the results to a file 
format, such as a text delimited file that is readable by 
5 a third-party data analysis tool. 

The system can include a use-case to define 
optimization constraints. To perform any simulation that 
requires a LP problem to be solved, the user must specify 
the constraints (upper and lower bounds) placed on all 

10 the reactions in the network and provide an objective 

function. These constraints define the conditions that 
are being simulated, such as growth phenotype under 
aerobic or anaerobic conditions or with glucose or 
without glucose. This use-case interacts with the user by 

15 providing a list of reactions and associated constraints 
from which a user can view and modify constraint values. 
Often times there are common constraint sets that the 
user will continuously use. So as not to require the 
user to repetitively enter common constraint sets, the 

20 system can store predefined constraint sets for 

particular models that are defined as templates from 
which a user can select and load one that is desired. 
Thus, this use-case provides a user with the option to 
select and load a predefined constraint template or 

25 select a constraint set from a previous simulation to use 
as the starting conditions, which can then be modified 
and used immediately or saved for future use. 

The system can include use-cases to perform any 
of a number of optimizations. A use-case to perform a 
30 single optimization is used to calculate a single flux 
distribution demonstrating how metabolic resources are 
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routed as determined from the solution to one LP problem. 
A use-case to perform a deletion analysis is used to 
calculate the consequences of deleting at least one gene, 
at least one biomolecule, or at least one reaction and 
5 running multiple LPs for each deletion case. A use-case 
to perform a robustness analysis is used to assess the 
effects of reducing the allowed activity through a 
particular metabolic reaction leading to a series of LP 
problems solved at each of the activity levels within a 

10 range. A use-case to perform a phase plane analysis is 
used to calculate the range of characteristic functions 
that a network can display as a function of variations in 
the activity of multiple reactions wherein an LP problem 
is solved for every combination of parameters. A use- 

15 case to perform a time-course analysis is used to analyze 
the transient shifts that occur in a network over a time 
period wherein an LP problem is solved at each time 
point . 

The use-cases for the various simulation types 
20 include features that allow access to linear programming 
algorithms and selection of parameters and data to be 
analyzed by the linear programming algorithms. These 
features include, for example, menus to load a network 
model, set constraints on all reactions and select an 
25 objective function. A simulation type use-case can have 
a user interface . that includes a main series of panels 
containing all of the intra-system reactions, input or 
output exchange reactions, demand exchange reactions, and 
temporary reactions that have been selected for a 
30 particular simulation. Upper and lower bound constraints 
on reactions can be specified by a user, for example, by 
changing the constraints displayed in a panel on the user 
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interface. Additionally , the user can select any 
reaction to be set as an objective function (such as a 
reaction representing cellular growth, ATP production, or 
a particular enzymatic reaction) . 

5 Results from each of the simulations can be 

viewed by a use-case of the simulation module. 
This use-case enables the user to view result data for a 
single optimization. Once a simulation has been run the 
solution can be output to a graphical user interface in 

10 any of a variety of acceptable formats for displaying 

simulation results including, for example, a table format 
or on a map. For any linear programming problem there 
are two sets of solutions, the primal solution and dual 
solution. Both the primal solution consisting of the 

15 flux values of all the reactions and the dual solution 
containing the reduced costs for the reactions and the 
shadow prices of the metabolites can be displayed. 

A use-case for comparing simulation results is 
also provided and can be used to simultaneously view 
20 tables or graphs from multiple simulations. A use-case 
is also provided for comparing simulation results to 
empirical results using similar tabular or graphical 
outputs . 

A robustness analysis can be performed by 
25 selecting a particular reaction or set of reactions for 
which the allowable flux level is reduced and running a 
simulation with the flux for the reaction (s) reduced 
using the use-case for performing a robustness analysis. 
From this use-case a user can select one or more 
30 reactions and then specify a set of constraints on the 
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reaction (s) or, in the case where incremental changes in 
constraints are to be analyzed, a step size increment by 
which the constraints will be changed can be set. The 
results of the simulation can be output to the graphical 
5 user interface in a tabular or graphical form using the 
use-case for viewing results of a robustness analysis. 

A phase plane analysis can be performed by 
calculating phase planes based on user defined parameters 
for particular reaction variables and value ranges. Here 

10 again the user specifies underlying constraint conditions 
and an objective value from the use-case for performing 
the simulation. The system runs all of the required 
single optimizations for one simulation and the results 
are presented using the viewing use-case in, for example, 

15 a tabular format or in a graphical representation. 

Following the simulations a shadow price analysis is 
performed to identify the different phases within the 
parameter space along with the isoclines for particular 
reactions specified by the user. As in all of the 

20 simulation type use-cases a particular point (or single 
optimization) can be selected and the system will 
generate the detailed solution of the corresponding 
single optimization for further analysis. 

Another simulation type is the time-course 
25 analysis which is performed to simulate transient 

cellular responses. In the use-case . for performing time- 
course analysis the user selects the baseline constraints 
and initial conditions from which to begin the 
simulation. The changes in extracellular reactant 
30 concentrations are calculated as a function of the 
uptake/secretion rate of the reactant, an initial 
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concentration, and the time increment specified by the 
user. The results can be viewed in a table or on graph 
charting the changes to the parameters in the analysis as 
a function of time using the use-case for viewing time- 
5 course analysis results. 

A data mining module can be included which 
provides the ability to evaluate the content of the 
models that have been developed. A wealth of knowledge 
can be derived from simple queries of the model content 

10 that do not necessarily rely on the simulation 

capability. A data mining module is available to manage 
all of these non-simulation related analyses. This 
includes the ability to ask questions concerning the 
reactions , proteins , and genes in various models. The 

15 • focus can be placed on one model in particular or on 

comparisons between many models. Text based or map based 
comparisons and result analysis are available. 
Metabolite connectivity studies can be performed as well 

A data mining module provides a number of use- 
20 cases to view data stored in various data bases, models 
or results files. A use-case to view an atlas allows a 
user to study network models by browsing through a set o 
network diagrams or maps. Similarly, a use-case to view 
model content allows a user to evaluate the content of 
25 the models using features such as browsing gene, protein 
and reaction related information in a tabular form, 
viewing model content on reaction maps or viewing 
gene-protein-reaction associations in a graphical 
association diagram. A reaction data base or compound 
30 database can be evaluated using use-case to view each. 
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A use-case can be included to perform a general 
content search of models. It includes the ability to ask 
questions concerning the reactions, proteins, and genes 
with the option to search within one model or across all 
5 models. Models can also be evaluated using a use-case to 
compare model content which allows a user to produce 
comparisons between many models using text-based or map- 
based comparisons and result analysis. 

Connectivity of reactants in a model can be 
10 evaluated using a use-case provided by the invention. 
This use-case includes the ability to view reactant 
occurrences on a map, view the connectivity for a 
particular reactant or a model in a tabular form or in 
terms of a connectivity graph. 

15 The genetic content of a bioparticle can be 

viewed using a use-case of the invention. This use-case 
includes features such as the ability to browse a gene 
index, view basic genetic content or view 
gene-protein-reaction associations . 

2 0 A number of additional modules also can be 

included in a computer system of the invention. These 
modules include, for example, an experimental data module 
for the integration and analysis of experimental data 
sets from high throughput experimental technologies such 

25 as gene expression arrays, protein expression arrays, 
protein-protein interaction arrays or metabolite 
profiling. Within this module experimental data sets can 
be compared against simulation results and enable the 
user to take advantage of experimental information for 

30 the iterative improvement of the model content and its 
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predictive capabilities. In addition to the experimental 
data module a gene sequence analysis module can be used 
to manage the process of annotating genomes to generated 
updated gene indices that are used to support model 
construction efforts. A pathway design module can also 
be introduced to allow for the network models to meet 
certain production requirements that a metabolic engineer 
may be seeking to design in a bacteria. This module also 
allows for the calculation of extreme pathways and 
related types of calculations which focus on the 
structural aspects of the metabolic networks that make up 
individual in silico models. 



EXAMPLE I 

Associating Genes, Proteins , and Reactions 

This example describes construction of a 
network model and a reaction index for the network model. 
This example demonstrates interactions of a user with the 
model content browser to associate the chosen ORFs to 
protein, and proteins to reactions. This example further 
demonstrates how this information is modeled from an 
object perspective and a data schema. 

A reaction index was constructed to include 
reaction components for both gene-associated and non 
gene-associated reactions . Gene-associated reactions 
were added to the reaction index as follows. 
Associations in the reaction index were formed based on 
known or putative associations of a reaction to the 
proteins or enzymes which enable or catalyze the reaction 
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and the open reading frames (ORFs) that code for these 
proteins. The associations were formed to capture the 
relationship between the reactions and proteins as well 
as between the proteins and ORFs such that connectivity 
5 between the reaction, protein (s) enabling the reaction 
and ORF(s) encoding the protein. 

The associations formed in the reaction index 
were displayed for review and evaluation by a user. The 
first panel of Figure 7 shows a display of the 

10 association in which one ORF (b2779) encodes one protein 
(Eno) which catalyzes one reaction (ENO) . Non-linear 
associations were also formed and displayed so as to 
capture the logic within the association. The non-linear 
associations for the PYRDH reaction are shown in the 

15 second panel of Figure 7 , where the requirement for both 
the b0114 and b0115 ORFs to encode the AceEF protein is 
indicated by the "AND" logic operator. Another non- 
linear association that was formed and displayed was that 
shown in the third panel of Figure 7 where two proteins 

20 (Tkt-1 and Tkt-2) encoded by separate genes (b2935 and 
b2465, respectively) are each capable of enabling the 
same two reactions (TKT1 and TKT2). The fourth panel of 
Figure 7 shows a display of the associations formed for 
the G3PDH reaction can be catalyzed by either the GapC or 

25 GapA protein, the former being encoded by two ORFs (bl416 
and bl417) and the latter being encoded by a single ORF 
(bl77 9) . The "OR" relationship between the GapC and GapA 
isozymes is displayed by multiple lines to the same 
reaction. 

30 The displays shown in Figure 7, by modeling 

associations, allowed evaluation of the network model and 
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its constituent reaction components at the gene, protein, 
or reaction level or at a combination of all three 
levels- In constructing the network model the 
associations were evaluated to determine the effects of 
5 adding or eliminating a reaction component at one level 
upon reaction components at another level. By viewing 
the associations shown in the third panel of Figure 7 , it 
was determined that removal of either the b2935 or b24 65 
ORF from the network model did not prevent flux through 

10 the TKT1 or TKT2 reactions. The association diagram 

displayed in the fourth panel of Figure 7 indicated that 
presence of either the bl779 ORF or the combination of 
the bl416 and bl417 ORFs will allow flux to occur through 
the G3PDH reaction. Thus, changes at the genetic level 

15 were readily correlated to biochemical activity of 
associated proteins and their reactions. 

In the course of forming associations, for each 
reaction, the identity of proteins required or capable of - 
performing the reaction was determined. For each 

20 protein, the number of subunits required for activity of 
the protein was determined. For each subunit, the number 
of ORFs that encode the subunit was determined. During 
iterative model construction, associations were formed 
and based upon display of the associations reaction 

25 components were evaluated for inclusion in the model. 

The gene-protein-reaction associations were 
formed in the Model Content Browser during the course of 
constructing the in silico network model. The Model 
Content Browser was accessed from the Model Construction 
30 main window by selecting the "Model Content Browser" 

button from the vertical toolbar shown in Figure 8. The 
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system opened the Model Content Browser window and 
displayed the gene index for the organism linked to the 
loaded model edition. 

The process of adding a gene-associated 
5 reaction to a model was divided into the following two 
steps. First, ORF-protein associations were formed. 
Second, protein-reaction Associations were formed. in 
the first step, one or more ORFs that should be 
associated with a reaction were identified. The gene 

10 index for the bioparticle was displayed as shown in 

Figure 9. The user navigated through the index using the 
slider bars that flank the index display. Once 
identified an appropriate gene was selected by activating 
the option "include" from a pop-up menu, as shown for the 

15 b0114 and b0115 ORFs in Figure 9. The selected ORFs were 
automatically added to the GENE-Proteln-Association 
Properties panel shown in the upper right portion of the 
screen shown in Figure 9. 

20 After selecting the b0114 and b0115 ORFs, an 

association was formed with the protein they encode. As 
shown in the upper right portion of the screen shown in 
Figure 10 , the AceEF protein was entered into the 
"Protein" entry field, thereby being associated to the 

25 b0114 and b0115 ORFs. The protein was selected from a 

drop-down list for the "Protein" entry field. If desired 
the protein' s abbreviation can be manually typed into the 
entry field. The system sent an automatic query to 
determine if the protein already existed in the system. 

30 Because the AceEF protein did exist the protein's name 
was populated in the field below the "Protein" entry 
field (see Figure 10) . In cases where the protein does 
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not exist, then the system enables an entry field where 
the user can enter the protein's full name. 

Once the ORF-protein association was correctly 
entered into the appropriate fields by the user, the 
5 apply button was clicked, in order to form the ORF- 
Protein association in the network model. The system 
responded by creating the appropriate database records 
and displayed the created associations visually in a 
graphical association viewer as shown in the lower right 
10 corner of the screen of Figure 11. 

The information describing the association was 
stored in a series of relational database tables. The 
following database records were created for the (b0114 
and b0115) — AceEF association of Figure 11. A peptide 
record was created containing the amino acid sequence of 
the polypeptide. In this case, the amino acid sequence 
was translated from the b0114 and b0115 ORFs . The 
peptide record was linked to the gene records for the 
aceE and aceF ORFs. Also created was a 

PeptideProteinAssociation record which represented the 
"AND" association of ORFs "b0114" and "b0115" to protein 
"AceEF". Further two PepPepProteinAssociation records 
were created to link ORFs "b0114" and "b0115" to the 
"AND" association record. These records entered as set 
forth above with respect to Figure 11 was stored in the 
proper database according to the object model shown in 
Figures 3 and 4. 

As set forth above in relation to Figure 11, 
multiple genes had to be associated with one protein in 
30 an "AND" relationship. The "AND" relationship was 



15 
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established automatically by the system upon the user 
entering the relationship in the "Gene Protein 
Association Properties" panel and sending the "apply" 
command. As shown in Figures 11 and 12, the graphical 
5 viewer represents this type of association with an 

symbol. An AND relationship between multiple genes and a 
protein reflects the quaternary structure of the protein 
including multiple subunits. 

There are two isozymes of the AceEF protein 
both capable of performing the PYRDH reaction. The first 
isozyme is encoded by the b0114 and b0115 genes. The 
second isozyme is encoded by the b2095 ORF. The 
relationship of the isozymes to the reaction was captured 
with an "OR" logic operator. As shown in Figure 13, the 
graphical association viewer represents an "OR" 
association by drawing multiple lines between the ORFs 
and the protein. The "OR" association is established 
when the user associates ORFs separately with the same 
protein. 

Next associations were formed between proteins 
and reactions. The Protein Index view was accessed by 
clicking on the "Protein Index" tab in the Model Content 
Browser. The system displayed all proteins that are 
associated with the model in a table as shown in Figure 
25 14. The appropriate protein, in this case AceEF, was 

selected from the protein index via the "Include" option 
from a pop-up menu as shown in Figure 15. In response 
the system populated the selected protein in the Protein- 
Reaction Association Properties panel on the right side 
30 of the screen. 



10 
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A reaction associated with the aceEF protein 
was entered into the "reaction" field. In this case the 
system found the reaction based on the abbreviation 
entered and populated the full name and equation in the 
5 appropriate fields. If the user does not know the 

reaction's abbreviation, the x \./' button can be selected 
to open a Reaction Browser window where reactions can be 
looked up from the reaction database based on any of a 
number of various criteria. Once the association was 

10 correctly entered the "apply" button was clicked to form 
the Protein-Reaction association in the network model. In 
response, the system then created the appropriate 
database records and displayed the created association ( s ) 
visually in a graphical association viewer located in the 

15 lower right corner of the screen shown in Figure 16. 

The system created the following database 
records for associations formed as described above in 
relation to Figure 16. A ModelReaction record was 
created to link the chemical reaction to the model. A 
20 ProteinReactionAssociation record was created to link 
the protein "AceEF" to the model reaction. A 
ProtProtReactionAssociation record was created to link 
the ProteinReactionAssociation to protein "AceEF". 

Protein-reaction "AND" and "OR" associations 
25 were established and displayed essentially as set forth 
above in regard to ORF-protein associations. A display 
of a protein-reaction "AND" association is shown in the 
graphical viewer in the lower right hand corner of the 
screen shown in Figure 17. A display of a protein- 
30 reaction "OR" association is shown in the graphical 
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viewer in the lower right hand corner of the screen shown 
in Figure 18 . 



describing a particular reaction are available and have 
5 been entered into the reference database, a link is 

provided to the reference by a "book icon" in the left 
hand column. For the reaction list shown on the display 
of Figure 17, the ACTL, AKGDH and PCK reactions have 
links to references. 

10 Figure 17 also shows a display in which the 

model reaction properties viewer is opened. In this 
viewer is shown information related to the confidence 
rating of the selected reaction. An overall confidence 
score is provided as well as a table showing the 

15 confidence details for five different categories. 

Confidence details and confidence scores are described in 
Example II. 

As shown in Figures 3 and 4, the following 
classes participate in the creation of Gene-Protein 
20 Associations: 



As shown in Figure 17, where references 



(1) 
(2) 
(3) 



Peptide 

Peptide ProteinAssociat ion and 
Protein . 
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The following classes participate in the creation of 
Protein-Reaction Associations: 



(1) Protein 

(2) ProteinReactionAssociation and 
5 (3) ModelReaction. 



The following tables participate in the creation of Gene- 
Protein Associations: 

(1) Peptide 

(2 ) PeptideProteinAssociation 

10 (3) PepPepProteinAssociation and 

(4) Protein. 



The following tables participate in the 
creation of Protein-Reaction Associations: 



(1) Protein 

15 (2) ProteinReactionAssociation 

(3) ProtProtReactionAssociation and 

(4 ) ModelReaction. 



EXAMPLE II 

Heuristic Algorithm for Confidence Level 



20 This example demonstrates a heuristic algorithm 

for determining overall confidence for inclusion of a 
reaction component in a particular network model based 
upon the level of information acquired in each of five 
categories . 
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The confidence levels range on a scale from 
zero to four with four being the highest rating level. A 
simple five level scale is adequate to distinguish 
between reactions with low confidence versus those with 
5 high confidence. The algorithm takes the level of 

significance assigned to each information category and 
filters them into a quantitative confidence level. The 
five levels will provide a basic indication of the 
confidence that the model content developer has in a 
10 reaction and the associated protein (s) and ORF(s) being 
included in a model. The meaning of the five levels is 
provided below. 

Level 0 - the reaction has no calculated 
confidence 

15 Level 1 - the reaction is supported by minimal 

evidence or even no evidence 

Level 2 - the reaction is supported by a fair 
amount of evidence 

Level 3 - the reaction is highly probable with 
2 0 ample evidence 

Level 4 - the reaction is certain to occur and 
has been validated 

The algorithm is based on the following equation: 

5 

CV = X InfoType i * InfoLevel f 

25 where CV is the confidence value that will be used to 
determine the confidence level, InfoTypei refers to a 
preset numerical value established for each of the five 
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information types (biochemical, genetic, genomic, 
physiological, modeling) , and InfoLevel ± refers to the 
preset numerical value for the information level that is 
associated with the specific information type. 

5 The following values were used for the preset 

numerical values for the information type and level: 

InfoType 

Biochemical 10 

Genetic 8 

10 Genomic 5 

Physiological 3 

Modeling 1 

Inf olevel 

Not evaluated 0 

15 None 0.1 

Low 1 

Medium 2 

High 3 
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Table 1 provides the range of confidence values 
that will correspond to the confidence levels to be 
prescribed to each of the reactions. 



Table 1 



Confidence 


Value Range 


Confidence 


Lower Value 


Upper Value 


Level 


0 


0 


0 


0.1 


8 


1 


8.1 


16 


2 


16.1 


22 


3 


22.1 


81 


4 



This framework for calculating the confidence 
rating allows for future alterations to the preset 
numerical values and ranges associated with each of the 
15 different information levels and types based on 

experiences gathered from implementing the confidence 
rating system described above. 

EXAMPLE III 

Identification and satisfaction of a macro requirement 
2 0 deficiency 

This example describes analysis of a network 
model to identify a gap. 



The user selects the "Run Gap Analysis" button 
to activate the network analysis. In response, the 
25 system activates the network analysis and identifies the 
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presence of gaps defined as either metabolites that occur 
only once as a reactant or metabolites that occur 
multiple times as only a substrate or product with all 
the reactions being irreversible. 

These situations will cause the associated 
reactions never to be utilized in the model simulations. 
For each gap, the system displays the name of the 
compound, the compartment in which the compound occurs, a 
description that indicates if the compound is consumed 
only or produced only, a checkbox that enables users to 
indicate which gaps have been reviewed. All gaps are 
sorted by compound abbreviation. 

In the following two examples, A and B occur 
only once as a reactant. A and B represent gaps if the 
15 reaction is reversible or irreversible. 

B 

A<^ B 

In the following example, B occurs multiple 
times as only a product (B is produced only) and all 
20 reactions it participates in are irreversible. B 
represents a gap. 

A~> B<r~ C 

In the following example, B occurs multiple 
times as only a substrate (B is consumed only) and all 



5 



10 
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reactions it participates in are irreversible. B 
represents a gap. 

A<r- B-+ C 

In the following example, B occurs multiple 
5 times as only a product (assuming that the second 

reaction was expressed as C <-> B and not as B <-> C) but 
one reaction is reversible. B does not represent a gap. 

C 

Throughout this application various 
10 publications have been referenced within parentheses. 

The disclosures of these publications in their entireties 
are hereby incorporated by reference in this application 
in order to more fully describe the state of the art to 
which this invention pertains. 

15 The term ^comprising" is intended herein to be 

open-ended, including not only the recited elements, but 
further encompassing any additional elements. 

Although the invention has been described with 
reference to the disclosed embodiments, those skilled in 

20 the art will readily appreciate that the specific 
experiments detailed are only illustrative of the 
invention. It should be understood that various 
modifications can be made without departing from the 
spirit of the invention. Accordingly, the invention is 

25 limited only by the following claims. 
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What is claimed is: 

1. A computer implemented process for 
constructing a scalable output network model of a 
bioparticle, comprising the computer implemented steps 

5 of: 

(a) accessing a database of network gene 
components comprising an annotated network set of open 
reading frames (ORFs)of a bioparticle genome; 

(b) forming a data structure associating said 
10 network gene components with network reaction components, 

said data structure establishing a data set specifying a 
network model of connectivity and flow of said network 
reaction components, and 

(c) transforming said data set into a 

15 mathematical description of reactant fluxes defining said 
network model of connectivity and flow, wherein said 
mathematical description defines a scalable output 
network model of a bioparticle. 

2. The process of claim 1, wherein forming 
20 said data structure further comprises: 

(a) selecting an ORF from said annotated 
network set encoding a gene product having a network 
reaction function; 

(b) forming a data structure comprising said 
25 selected gene product, said data structure associating 

network gene components and network reaction components 
comprising cognate ORFs, encoded gene products, network 
reactions and reaction constituents, and 

(c) repeating steps (a) and (b) selecting 
30 another ORF from said annotated network set until 

substantially all of said network gene components of said 
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annotated network set have been surveyed for encoding a 
gene product having a network reaction function to 
produce a data structure establishing a data set 
specifying a network model of connectivity and flow. 

5 3. The process of claim 2, further 

comprising: 

(a) determining the occurrence of a 
constituent gene product for said selected encoded gene 
product ; 

1° (b) determining the occurrence of an 

additional gene product participating in said network 
reaction; 

(c) determining the occurrence of an 
alternative network reaction exhibited by a surveyed gene 

15 product; 

(d) incorporating identified constituent gene 
products, participating gene products or alternative 
network reaction into said data structure. 

4. The process of claim 1, further comprising 
20 incorporating a network reaction that is not gene-encoded 

and corresponding reaction constituents into said data 
structure . 

5. The process of claim 1, further comprising 
elemental balancing on at least one network reaction. 

25 6- The process of claim 1, further comprising 

charge balancing on at least one network reaction. 
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7. The process of claim 1, further comprising 
incorporating an exchange reaction for an external 
reaction component and corresponding reaction 
constituents into said data structure. 

5 8. The process of claim 7, wherein said 

external reaction component comprises a metabolite or a 
biochemical demand constituent. 

9. The process of claim 8, wherein said 
biochemical demand further comprises an aggregate 

10 reactant demand flux defining a phenotypic output for 
growth . 

10. The process of claim 9, wherein said 
phenotypic output for growth comprises biomass 
production . 



15 11. The process of claim 8, wherein said 

biochemical demand further comprises an aggregate 
reactant- demand flux defining a- phenotypic output . 
selected from the group consisting of energy production, 
redox equivalent production, catabolite production, 

20 biomass precursors, polypeptide production, amino acid 
production, purine production, pyrimidine production, 
lipid production, fatty acid production, cofactor 
production, production of a cell wall component and 
transport of a metabolite. 



25 12. The process of claim 1, wherein said data 

structure comprises reactants, products and 
stoichiometric coefficients . 
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13. The process of claim 1, wherein said 
mathematical description comprises linear equations and 
inequalities . 

14. The process of claim 13, wherein said 
5 mathematical description comprises a stoichiometric 

matrix. 

15. The process of claim 13, wherein said 
mathematical description comprises differential 
equations . 

10 16. The process of claim 1, further comprising 

calculating a phenotypic output of said network model 
from said mathematical description. 

17. A computer implemented process for 
constructing a scalable phenotypic output network model, 
15 comprising the computer implemented steps of: 

(a) accessing a database of network gene 
components comprising an annotated network set of open 
reading frames (ORFs)of a bioparticle genome; 

(b) forming a data structure associating said 
20 network gene components with network reaction components, 

said data structure establishing a data set specifying a 
network model of connectivity and flow of said network 
reaction components; 

(c) modifying said data set to enumerate a 
25 biochemical demand on said specified network model, and 

(d) transforming said modified data set into a 
mathematical description of reactant fluxes defining said 
network model of connectivity and flow, wherein said 
enumerated biochemical demand corresponds to an aggregate 
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reactant demand flux defining a phenotypic output of said 
network model of a bioparticle. 

18. The process of claim 17, wherein forming 
said data structure further comprises: 
5 (a) selecting an ORF from said annotated 

network set encoding a gene product having a network 
reaction function; 

(b) forming a data structure comprising said 
selected gene product, said data structure associating 

0 network gene components and network reaction components 
comprising cognate ORFs, encoded gene products, network 
reactions and reaction constituents, and 

(c) repeating steps (a) and (b) selecting 
another ORF from said annotated network set until 

5 substantially all of said network gene components of said 
annotated network set have been surveyed for encoding a 
gene product having a network reaction function to 
produce a data structure establishing a data set 
specifying a network model of connectivity and flow. 

20 19. The process of claim 18, further 

comprising : 

(a) determining the occurrence of a 
constituent gene product for said selected encoded gene 
product; 

25 (b) determining the occurrence of an 

additional gene product participating in said network 
reaction; 

(c) determining the occurrence of an 
alternative network reaction exhibited by a surveyed gene 
30 product, and 
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(d) incorporating identified constituent gene 
products, participating gene products or alternative 
network reaction into said data structure. 

20. The process of claim 17, further 

5 comprising incorporating a network reaction that is not 
gene-encoded and corresponding reaction constituents into 
said data structure. 

21. The process of claim 17, further 
comprising elemental balancing on at least one network 

10 reaction. 

22. The process of claim 17, further 
comprising charge balancing on at least one network 
reaction . 

23. The process of claim 17, further 

15 comprising incorporating an exchange reaction for an 

external reaction component and corresponding reaction 
constituents into said data structure. 

24. The process of claim 23, wherein said 
external reaction component comprises a metabolite or a 

20 biochemical demand constituent. 

25. The process of claim 17, wherein said 
biochemical demand further comprises an aggregate 
reactant demand flux defining a phenotypic output. 

26. The process of claim 25, wherein said 

25 phenotypic output further comprises an aggregate reactant 
demand flux defining growth. 
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27. The process of claim 25 , wherein said 
phenotypic output further comprises biomass production. 

28. The process of claim 17, wherein said 
biochemical demand further comprises an aggregate 

5 reactant demand flux defining a phenotypic output 

selected from the group consisting of energy production, 
redox equivalent production, catabolite production, 
biomass- -precursors, polypeptide production, amino acid 
production, purine production, pyrimidine production, 
10 lipid production, fatty acid production, cofactor 

production, production of a cell wall component and 
transport of a metabolite. 

29. The process of claim 17, wherein said data 
structure comprises reactants, products and 

15 stoichiometric coefficients. 

30. The process of claim 17, wherein said 
mathematical description comprises linear equations and 
inequalities . 

31. The process of claim 30, wherein said 
20 mathematical description comprises a stoichiometric 

matrix . 

32. The process of claim 30, wherein said 
mathematical description comprises differential 
equations . 

25 33. The process of claim 17, further 

comprising calculating a phenotypic output of said 
network model from said mathematical description. 
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34. A computer implemented process for self- 
optimizing a network model of a bioparticle, comprising 
the computer implemented steps: 

(a) accessing a database of network gene 

5 components comprising an annotated network set of open 
reading frames (ORFs)of a bioparticle genome; 

(b) forming a data structure associating said 
network gene components with network reaction components, 
said data structure establishing a data set specifying a 

10 network model of connectivity and flow of said network 
reaction components ; 

(c) transforming said data set into a 
mathematical description of reactant fluxes defining said 
network model of connectivity and flow; 

15 (d) determining the competence of said 

connectivity and flow within said network model, said 
competence indicating underinclusion or overinclusion of 
network reaction component content of said network model, 
and 

20 (e) identifying an ameliorating network 

reaction component capable of augmenting said competence 
of said network model, incorporation of said ameliorating 
network reaction component into said data structure 
producing a modified data structure specifying in an 

25 optimized network model of said bioparticle. 

35. The process of claim 34, wherein said 
network comprises a metabolic network. 



30 



36. The process of claim 35, wherein said 
metabolic network further comprises a plurality of 
network pathways of a bioparticle genome. 
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37. The process of claim 34 , wherein forming 
said data structure further comprises: 

(a) selecting an ORF from said annotated 
network set encoding a gene product having a network 

5 reaction function; 

(b) forming a data structure comprising said 
selected gene product, said data structure associating 
network gene components and network reaction components 
comprising cognate ORFs, encoded gene products, network 

10 reactions and reaction constituents, and 

(c) repeating steps (a) and (b) selecting 
another ORF from said annotated network set until 
substantially all of said network gene components of said 
annotated network set have been surveyed for encoding a 

15 gene product having a network reaction function to 
produce a data structure establishing a data set 
specifying a network model of connectivity and flow. 

38. The process of claim 37, further 
comprising: 

20 (a) determining the occurrence of a 

constituent gene product for said selected encoded gene 
product; 

(b) determining the occurrence of an 
additional gene product participating in said network 

25 reaction; 

(c) determining the occurrence of an 
alternative network reaction exhibited by a surveyed gene 
product, and 

(d) incorporating identified constituent gene 
30 products, participating gene products or alternative 

network reaction into said data structure. 
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39. The process of claim 34 , further 
comprising incorporating a network reaction that is not 
gene-encoded and corresponding reaction constituents into 
said data structure. 

5 40. The process of claim 34, further 

comprising elemental balancing on at least one network 
reaction . 

41. The process of claim 34, further 
comprising charge balancing on at least one network 

10 reaction. 

42. The process of claim 34, further 
comprising incorporating an exchange reaction for an 
external reaction component and corresponding reaction 
constituents into said data structure. 

15 43. The process of claim 34, further 

comprising incorporating a biochemical demand into said 
data structure. 

44. The process of claim 34, further 
comprising : 

20 (a) determining the occurrence of a network 

reaction component satisfying a macro requirement 
deficiency in structural architecture of said network 
model, and 

(b) incorporating an identified network 
25 reaction component satisfying said macro requirement 

deficiency into said data structure to supplement said 
connectivity and flow of said network model. 
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45. The process of claim 34, further 
comprising executing a heuristic logic decision algorithm 
determining confidence of said network reaction 
components within said data structure. 

5 46. The process of claim 34, wherein said 

mathematical description comprises linear equations and 
inequalities. 

47. The process of claim 46, wherein said 
mathematical description comprises a stoichiometric 

10 matrix. 

48. The process of claim 46, wherein said 
mathematical description comprises differential 
equations . 

49. The process of claim 34, further 

15 comprising determining said competence by solving said 
mathematical description for a single optimization 
solution, wherein the ability of said network model to 
produce a pathway flux indicates a competent network 
reaction component content. 



20 50. The process of claim 49, further 

comprising solving said mathematical description for a 
plurality of single optimization solutions. 
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51. A computer implemented process for 
constructing a data structure specifying a network model 
of a bioparticle, comprising the computer implemented 
steps : 

5 (a) accessing a database of network gene 

components comprising an annotated network set of open 
reading frames (ORFs)of a bioparticle genome; 

(b) selecting an ORF from said annotated 
network set encoding, a gene product having a network 

10 reaction function; 

(c) determining the occurrence of a 
constituent gene product for said selected encoded gene 
product; 

(d) determining the occurrence of an 

15 additional gene product participating in said network 
reaction; 

(e) forming a data structure from said 
selected and determined gene products , said data 
structure associating said network gene components and 

20 network reaction components comprising cognate ORFs, 
encoded gene products, network- reactions and reaction 
constituents, and 

(f) repeating steps (a) -(e) selecting another 
ORF from said annotated network set until substantially 

25 all of said network -gene components of said annotated 
network set have been surveyed for encoding a gene 
product having a network reaction function to produce a 
data structure establishing a data set specifying a 
network model of connectivity and flow. 
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52. The process of claim 51, further 
comprising the steps of: 

(a) determining the occurrence of an 
alternative network reaction exhibited by a surveyed gene 

5 product, and 

(b) incorporating an identified alternative 
network reaction and corresponding reaction constituents 
into said data structure. 

53. The process of claim 52, further 
10 comprising: 

(a) determining the occurrence of a 
constituent gene product or a gene product participating 
in said alternative network reaction, and 

(b) incorporating an identified constituent 
15 gene product or gene product participating in said 

alternative network reaction into said data structure. 

54. The process of claim 51, further 
comprising incorporating a network reaction that is not 
gene-encoded and corresponding reaction constituents into 

20 said data structure. 

55. The process of claim 51, further 
comprising elemental balancing on at least one network 
reaction . 



56. The process of claim 51, further 
25 comprising charge balancing on at least one network 
reaction. 
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57. The process of claim 51, further 
comprising incorporating an exchange reaction for an 
external reaction component and corresponding reaction 
constituents into said data structure. 

5 58. The process of claim 57, wherein said 

external reaction component comprises a metabolite or a 
biochemical demand constituent. 

59. The process of claim 51, further 
comprising incorporating a biochemicalx demand into said 

10 data structure. 

60. The process of claim 59, wherein said 
biochemical demand further comprises an aggregate 
reactant demand flux defining a phenotypic output of said 
network model. 

15 61. The process of claim 51, further 

comprising : 

(a) determining the occurrence of a network 
reaction component satisfying a macro requirement 
deficiency in structural architecture of said network 

2 0 model, and 

(b) incorporating an identified network 
reaction component satisfying said macro requirement 
deficiency into said data structure to supplement said 
connectivity and flow of said network model. 

25 62.. The process of claim 61, wherein said 

macro requirement deficiency comprises a pathway gap or a 
pathway dead-end. 



WO 03/106998 



PCT/US03/18838 



126 

63. The process of claim 62, further 
comprising identifying a singleton reactant. 

64. The process of claim 62, further 
comprising identifying a reactant participating solely in 

5 two or more irreversible network reactions. 

65. The process of claim 61, wherein said 
network reaction component comprises a substrate or a 
product . 

66. The process of claim 51, further 

10 comprising executing a heuristic logic decision algorithm 
determining confidence of said network reaction 
components within said data structure. 

67. The process of claim 66, wherein said 
inclusion of a network reaction component further 

15 comprises determining a confidence level from a 
hierarchical classification. 

68. The process of claim 67, wherein said 
hierarchical classifications are selected from the group 
consisting of biochemical, genetic, genomic, 

20 physiological and simulation modeling data. 

69. The process of claim 51, further 
comprising transforming said data set into a mathematical 
description of reactant fluxes defining said network 
model of connectivity and flow of network reaction 

25 components. 
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70. The process of claim 69, wherein said 
mathematical description comprises linear equations and 
inequalities . 

71. The process of claim 69, wherein said 
5 mathematical description comprises a stoichiometric 

matrix. 

72. The process of claim 69, wherein said 
mathematical description comprises differential 
equations . 

10 73. The process of claim 51, further 

comprising performing a validation test. 

74. A system for constructing a scalable 
output network model of a bioparticle, comprising: 

(a) an input data set of network gene 

15 components comprising an annotated network set of open 
reading frames (ORFs) of a bioparticle genome; 

(b) executable instructions forming a data 
structure associating said network gene components with 
network reaction components, said data structure 

20 establishing a data set specifying a network model of 
connectivity and flow of said network reaction 
components; 

(c) executable instructions determining the 
occurrence of a reaction component satisfying a macro 

25 requirement deficiency in structural architecture of said 
network model, inclusion of an identified reaction 
component satisfying said macro requirement deficiency in 
said data structure supplementing said connectivity and 
flow of said network model; 
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(d) a heuristic logic decision algorithm 
determining confidence of said network reaction 
components within said data structure, and 

(e) executable instructions mathematically 
5 describing from said data set reactant fluxes defining 

said network model of connectivity and flow, wherein said 
mathematical description defines a scalable output 
network model of a bioparticle. 



75. A system for constructing a scalable 
10 phenotypic output network model of a bioparticle, 
comprising: 

(a) an input data set of network gene 
components comprising an annotated network set of open 
reading frames (ORFs) of a bioparticle genome; 
15 (b) executable instructions forming a data 

structure associating said network gene components with 
network reaction components, said data structure 
establishing a data set specifying a network model of 
connectivity and flow of said network reaction 
2 0 components; 

(c) executable instructions modifying said 
data set to enumerate a biochemical demand on said 
specified network model, and 

(d) executable instructions mathematically 
25 describing from said modified data set reactant fluxes 

defining said network model of connectivity and flow, 
wherein said enumerated biochemical demand corresponds to 
an aggregate reactant demand flux defining a phenotypic 
output of said network model of said bioparticle. 
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76. A system for constructing a self- 
optimizing network model of a bioparticle, comprising: 

(a) an input data set of network gene 
components comprising an annotated network set of open 

5 reading frames (ORFs) of a bioparticle genome; 

(b) executable instructions forming a data 
structure associating said network gene components with 
network reaction components, said data structure 
establishing a data set specifying a network model of 

10 connectivity and flow of said network reaction 
components; 

(c) executable instructions mathematically 
describing from said data set reactant fluxes defining 
said network model of connectivity and flow; 

15 (d) executable instructions computing 

competence of said connectivity and flow within said 
network model, said competence indicating underinclusion 
or overinclusion of network reaction component content of 
said network model, and 

20 (e) executable instructions augmenting said 

competence of said connectivity and flow within said 
network model, said executable instructions specifying 
inclusion or exclusion of an ameliorating network 
reaction component, wherein incorporation of said 

25 ameliorating network reaction component into said data 
structure produces a modified data structure specifying 
an optimized network model of said bioparticle. 



WO 03/106998 



PCT/US03/18838 



1/18 



Client Oesktap 




Application Server 



Computational Engine 
Server 



:i M 



Database Server 




FIGURE 1 



WO 03/106998 



2/18 



PCT/US03/18838 



100 

Initiate Model Development 



no 



200 

Complete Model 

Requisition/ 
Approval Process 



300 
Collect All 
Relevant 
Organism Specific 
Information 




400 

Construct Network 
and Establish 
Model Content 
Associations 



500 
Preliminary 
Simulation Testing 
and Model Content 
Refinement 



600 

Review Content 
and Prelimiinary 
Simulation Results 
Analysis 




yes 



700 

Version Model and 
Document 



C 



Maintain Citation Library 



DBA/Bioinformatician 

p 
A 

Curator 

o 



Maintain Organism Specific Data A 

Model Developer 




o- 

Maintain Model Content 



Model Developer 



Validate Model 
Structure & Content Model Developer 



\ 



\ 



x <zz> 



Perform Model Testing A 

Model Developer 



o 



Validate Model 
Structure & Content 



A 

Model Developer 




800 

Model Versioning Complete 



Version & Release Model A , 

Model Developer 



FIGURE 2 



WO 03/106998 



PCT/US03/18838 




3/18 



FIGURE 3 



WO 03/106998 



PCT/US03/18838 



iff ifi 



5 



II 



ill Si llWM mm 

s^isiiiiiiilliillillilli 




-If 

illflillll 

5 tstt a 

eg e 



sip 

III 



f 



Ills 



fliilllfip g 



«-fffJiia 



5-C W££ 



lliii 



iifpiJiiii 

II ||||||||| 



H 



4/18 



WO 03/106998 



PCT/US03/18838 



5/18 



(^Initiate 



410 

Initiate Network Construction 




414 

Create/Load Model 
Edition 




413 
Create Open 
Edition 



The system will run a pre- 
simulation validation lest to 
determine if all of the 
requirements for the model 
are in place to begin full 
testing. Simulation not begin 
without running this 
validation check. 



416 
Load Model 
Content 



420 

Evaluate Gene 
Index & add Gene 
associated 
Reactions 




Evaluate Gene Index 



Model Developer 



430 






Add Known 






Biochemical 






Reactions Nan- 
Gene Associated 






T 






440 






Add I/O and 






Biomass 






Constituent 






Exchange Fluxes 







sintain Reaction Index 



Curator 



450 

Create Aggregate 
Biomass Demand 
Exchange Fluxes 



460 
Identify and 
Evaluate Network 
Gaps 



470 
Introduce 
Confidence Levels 
for All Model 
Reactions 




Maintain Exchange Reactionst 



i 



Model Developer 



Complete Model Reconstruction 



. T- 

struction . /\ 



Model Developer 



480 

Perform Validation 
of Model Structure 
and Content 




Validate Model 
Structure & Content Model Developer 



490 \ 
,End Pre-Testing Model Development/ 



FIGURE 5 



WO 03/106998 



PCT/US03/18838 



6/18 




Initiate Gene Index Evaluation 



Reveiw Product 
Description fo r 
next unclassified 
gene 



Exclude Reaction 

with Reason 
"Unkown Function" 



Exclude Gene with 
reason "Non- • 
Metabolic" 



Exclude Gene with 
reason "Unclear 
Annotation" 



Model Developer 




Identify Additional 
Genes for the 
subunits 




Add Protein to the 
Protein Index 



Associate Gene(s) 
With 
corresponding 
Protein 



Curator 



yes 




Developer 



Complete Gene Index Evaluation 



Jon^ 



FIGURE 6 



WO 03/106998 



PCT/US03/18838 



7/18 






FIGURE 7 



WO 03/106998 



PCT/US03/18838 



8/18 



tag 



'^i3«j£^&dwdthl*'«»» MS-«Omp Edition « CT.fi AW»|ii. j[«S*a7«0<& tt*MM) 



Stwetuie 



ValiaaKer* 



f'lopmrtitg 



«io4jt} 



Modet ConteJJt Browser . ; . 



- " - 



"Gene Index Pinton kuto; r; 



Coordinate! Narne j&sneStmfea 




Proiem km 


Protein 




| 180ibOG01 




thr oneton leader peptide I (none) 


22 








i 337jb0002 


*m 


aspartokmase i homoWine tie .. i (none) 


821 








II 2.601 


ib0GD3 




bomos etine Mrtase I (none) 


311 








} 3,734 jb0004 




sthreoflirw synthase I (none) 


425} 








1 5.234 


ibOOOS 




[off. hypo!he(ic«i protein ; ; (none) 


99 








I 6,459 bODOS 




orf. byporheSca! protein { (none) 


259 








I 7.863 ! b0007 




ifri sen br; ru .un-r rfproU n (i <rr: 


477' 






i 8,239 


<b0G08 




transalSoiase B limine) 


318| 






! 9,308 mum 




r 5 f f *' ! i { n 5 ^r ! m off 3 i fl , 


198! 






\ 1 0,494 IbOOIt) 


*-*>a# 


orf. hvpotnebeat protein ? (none) 


_9™". 






If ft .356 bO'ai 1 




Ptf-<:j -t- o>7dc«e2yci:. :<:• . (none) 


238 








| 10J25;bOD12 


htgA 


• c lirvoi auTsto rot ipn 3321 . unr 


197 








i 11.785 : : bD013 


ort.'h -0U1E :rsi m- <n ^ione^ 


L 135 








I 12,183 :b001 4 




^chaperons H o 7 u D n «. fc , 0 svnt (no n - , 










1 * 1 C l;001'6 




ilurpef n wiihDng^ heat shoe (none) 


377 








I 15,445 i'bODI 6 


WiJ 


lSl86hypcKn-t - f j 1 Qionef 


I 371 














[3186 an5Ts421 r>, ottieti r 'f n.ie) 


103 








i 16.880 


boon. 




Ooi orotfiri inir-oetcf wtii rn?inL. (none) 


70 








! "17.489 


biV.'O 




iki^'HanrruCor, pH ^t^^ idenl (no v; 










I 18,715 


b0020 






382 










b0021 


>r.33.i 


is pjDieiniitPe Itnone) 


163 








| 20,508 


bflC22 




(51 pcejein JtisA frone) 


92 








Ij 2i t 078 


b0023 




303 tibosorr.&i cubumtptoi?ir ( S .Hnone) 


aB 








L™31,t?.l 


bPM< 


t?0024 


.it oib b aip ola] 1 l^nDne> 


73 








|| 21.407 


b0025 " 




pulative tegulaiot wmO 










1 2091 


b()026 


iites 


1 I't^oi - if fjs svr > h Mase l(m»ne> 


<m 








1 25,207 


b0027 ffepA 


probpoprcilein stpnal paptidase { \rtv 1 s 


1G5 








| 25,826 


bOQ28 


u'^'Mf Fl Bz-t-pa !t } Dpept 3 foun^/ 


150 








1 26 A 277 


bQ029 


torntui i{ i nnnpnt . ;po if (uv , « 


'iV 








[ 27.283 


b0030 




orf, 1 n c-fi'isi Lit y-olblr tyone) 


305 


"_ 






S 26^374 


bOQ31 




an ii.rJi, ilc insiered i h v r 


274 








{[ 28.651' 


1.0032 


<carA 


carbamoyi-phospj ale syntriaias (not c> 


353 








( 30,817 






carrsrrjovJ-phospbar© Suisse J.. J (non»> 


j,07* 








] ' 34.195 


b0034 




tTanicripiionat regulator ofcai ojgi„|(nbrie) 


167 








| 35,382 


b003f 




possible synmei ?ofcDfactcirfo ,(none) 


204 








I' 38.270 


b003B 


*?a/0 


carnilJns racerniaVe 1 (none)™™ 


?9S 








J 37,839 


b0037 




probable trotonobeia , ineicamin...i^none) 


523 










b0038 


Sca/8 


icajainede^faiase' Krions) 


406 






f 40,386" 


b0039 




ore cuK 1 faruT'ife operon o<jdor. (none) 


38f 






< ' . . 








BM*,fo;A "" 







FIGURE 8 



WO 03/106998 



PCT/US03/18838 



9/18 



^Genetic Circuit Simulation - Modef Oevetopmem i 



F8e m View loots Help _ - ^ ^ ..„.,_ '.'.„„ '. . J. — > - - ' 



Mod*! 
Content .'j! 



1 Mo«i«r 



1 f.teJel Content Browser ' 



Gc ii Index " P«ite ,t Intte Ri Action Inat x 



S*H»S 1 All 
1 


11 , 


m ... 












5'Coordinate 


Name 


jOene Symbol 1 


. , Putative ID ; - jCelluiar RoUPtotem ten...] 


Protein 


115,724 


b0106 




pi.iaiiv* mil orJ memorane p<^t in invoh troiv'i 


401 1 * 


117,089 


b0107 


^OfB „ 


putative mtfepr< i nbra protein invotv utie) 1 


462 1 


{17,648 


b0108 


\ppiSD 


pxsltofn peptidase dependent protein i (none) I 


1471 ~ 


118,645 


b0109 




quinolinate phosphoribosyrtranafetase \{rmn&) j 


298;; 


118,733 


b0110 


UmpO 


reouiaies ampC Hoooe) . 


1£4 




119.201 


bom 




renuiatesampC 


(none) 




285 




121^51 


b0112 


\aroP 


urpiBMi amino ai id transport protein (m.,i&) 


458! 




122.092 


D0113 


\pdtiR 


transcriptional regulator toi r ? j di l t ton«) 


255? 




1 ziffl^ 
125,695' 

""127,812 


bQ114 
L'01 15 
b0116 




lipoarnids debydragenaTe Associati 






1 888 




IfpdA . 






475; 


131,260 


b0117 




putative membrane protein j ~p V rf„dP 






6181 


131. ii': ■ 


b0118 


lacnB 


aconitate hydrase B 1 I. • , . , : 




S68j 


134.340 


b0119 




•ortitootheli^ 5*E^™5,,, 




1371 


^13^82, 


b0120 




S-adenosytmetblonins decarboxylase ?<none) 


265! 


136.484 


b0121 




spermidine svntriase^ puirescihe amin... (none) 


289| 


137,040 
13? ( 0fi3 


b0122 




u ' <*t l i itscs prot m ?<none!> 


' 157 1 


b0123 


""Sac* 


irf, hy\ thsli ai pro) r ^l^5 ne ?-. 


5171 


141,225 


b0124 ktcd 


(jiucose deihvTttoyenase Unone) 




141^19 


b0125 


ftp? 


hvpoxanthtne phasphortbosvitransferase | 


(none) 




183| 


1 42.670 


b0125 




put-, ive carbonic anhdrsse (EC 4.2,1 1) f.none). 


221 | 


'14 2,773 


b01 27 




m \ ~ihi - A 1 1 -tutu ! <i - nr^inf f nt 3fatrat (nf< e 


306 


143,707 


b0l28 




Oif hyp tit;; i p"~t^<n 


j^mme) 




257:; 


144.577 


b0129 


iputatfve PTS enzytne ii 8 component 


(none) 




1471 


145.081 


|b0130 




orf, h>/pothetica! protein i (none) 


410 


146.8J4 


b0131 




aspartate 1-decatbo^lase 


(none) 




127| 


146.968 


b0132 










301 




148,795 


b0133 


|paDG 


pantothenate s^fnthetase 


(none) 


284j 


148.801 


b0134 


JoanS 


3-metnvt-2 oxobutanoate hydro^ymethyttr, 


(none) 


285 i 


160,953 


fa0l35 




pu! au'- luntHl.T 


(none) 


413 




i6i^9ja 


b0136 




putatK'& timbrfa? protein 




139 




152,231 


b0137 




nuijttv? fmt'fKV protein 


(none) 


202 




152,854 


}b0138 




nutsiut (.nifanahts^ uro'em 


(none) 


204 




: 155,426 


b0138 


mii 


probaPie outer membrane porin pioteini... 


(none) 


B5B 




155,201 


b0140 




probable piltn chaperons simliar to PapO 


JnoneJ 


247 




155,883 


b0141 




putative nmbriaWike protein 


^none) 


1S5 




157,732 


b0142 




7 r di t rydro-6-}vy'd t /••rn^hv^.^rm- pyrop 


(none) 


160 




153,093 


b0143 




,< ■ , - po me t 1 


(none) 




455 




180.112ib0i44 


ryacfB 




'■■> V s ; 




309 




: N* 













| Locus Name 


i PutatrvelD ! 


b0114 


pyruvate dehydrogenase {deeaitfoi^a,.r 


P.0116 


pvruvaTe der^dtogenase Cdih^rollppy„ J i 










Hc<s*s 4 AppV ft 1 



FIGURE 9 



WO 03/106998 



PCT/US03/18838 



10/18 



mmmmm^ 



file View Tools Help , ^ 



Gofftunl 
Brr»»e*f 

If'-; 

Strijetuje 



■IComem Browser ' ; , " ' < ^ ' , . ^ -^ *tJl 



Grnoji,^^ Pratom imtex Reaction litile* 



S* Coordinate! Mama f sene Syrribot 



Puiatiyejp 



Cellular Soi.JPtotem tan.., 



1 1 5,724 jbOl 06 



1 1 7,088 in01 07 



117.543100108 



..jm645XhD1D9,_ 



j putative integral membrane protein invotv.>: <none) 



118.73311)0110 



113.2311150111 



\ p re ! jpj n p e p'i ti a s s dependentjroteln iiii 



qu ribs n t< r mi phoribosyllran ft ras 



(none) 



(none) 



1 .21,551 I bpl 12 



j regulates ampC 



(none) _ 



j<none) 



3ce£, 



Tarornatft- amino acid transport protein j<hone) 



f iypj i alt* ^l™drog< r>nc/<.\ ei ama*\l'iz 
tfvrLV.ley-f 1 wogejiase (dlhydrclip lyitra 



&0115 




ite dehydrogenase fdecarb oxyla. 



FIGURE 10 



WO 03/106998 



PCT/US03/18838 



11/18 



S&Genetic Circuit Simulation ~ Model Development 



RJb pi* View Tools Help • ^ „ ; . r — : .- 



Mod*) : 
Content 



M«fcfc6fk 
0*P* 



Prop«rti*5 



^ fjto0ei Comem 8rcwsei ' ■ '- - : . : - 



Gene index ]jip|S^ffi^J^|tfon^5e^j 



js 4 Coordinate 


ttame 


£ene Symbol 


Putative ID 


iCeliuiarRot 


Protein l&rCj 


Protein j 


( 115,724 


b01 06 fftOJC 


putative integral membrane protein invotv. 


intone) 


401 






117,099 


b0107 


hofB 


putative iniearal membrane protein invotv.J (none) 


462 






117.549 


b0108 


ppdD 


orellDin peptidase dependent protein 


<<none) 


147 






118.645 


[b0109 


roidO 


quinatm rte p, orpnonU _,0tr-p-ji-M - 


<(none) 


298 






118,733 


b0110 


•ampD 


ireputatesjampC 


f(non^2I^, 


ie'i'l 






119,281 


[00111 


[amp£ 


regulates ampc 


jjOjone) 


385 






121.651 


b0112 


STOP 


aromatic amino acid transport protein 


Hnons) 


458 






122,082 
™" 123 017 


pom iu 


3ghR_ I 


transcnpliomil regulatorfor pyruvate deny, (none) t 
1 tHii'tit-' tiPli aronw;We fiier arbors e 'muie,) 


255 


'AceEF'" 








f Mu*«ie dp^voroqena.t (dlbYdrGHpoyttra 


, (none): : i 


. 621 AceEF 




127.912 b01 16 . 


\todA 


lipoamide det ydrod^nase OM.^Dri). corn 


^fjune) 


475 






131.260 


5b0117 




putative membrane protein 


\ (none) 


618 






131.615 


IbOiiO 


\acnB 


3 snitati i -! ( r,< B 


tjfnoneji 


866 






134,340ib0119 


\yacL 


orf, hypothetical protein 


^none) 


137 






135.592ib0120 




S-adeiiosyjmtith!Dii'r>Pdscarboxytase 

spem ine vrtlh t - puirescTn 3nni i 


[(none). _ 


265 






-j— - 




mei . 


|^one) 


289 






137,040 ib01 22 " 


vacC 


orf, hypothetical protein 


limine) 


157 






137.983 b01 23 


lyacK 


orf, hypotheiica) protem 


{ (none) 


517 






141,225 


jb0124 


tied 


alucose dehydrogenase 


(none) 


797 






141.419 


ib0125 


m 


hvpoxan n lephi •.phoribosYitransfsras 


i^none) 


183 






I 142 370 JU0 126 


\yadF 


putative i iiti,!fi!C anhtirsse lEO < 2 1 1)^ 




221 






142,773 


i 00127 


\y&(fG 


putative ATP'bi'nding component of a ban 


J (nonef 


309 






143.702 


ib0128 


jyscW 


orf, hypothetical protein 


{(none) 


257 






144,577 b0129 




putative PTS enzyme if B cornponent 


] (none) 


147 






I 145.081 1 bOI 30 


yatf£ 


orf hypt ttiejr *\ protein 


I (none) 


410 






j 146,694 


ib0131 


ipand 


[aspartate I'decarboxviase 


i (none) 


127 






146.968 


Ib0132 


{ yzdD 


on". h>pothetitai protein 


l(none) 


301 






143.795 jb0 133 


'•pane 




^none) 


284 






149,601 


5b0134 


ipanB 


3-n.HiTv.-; n'G^uS-.n ntp nvdit.i-anctti^tf 


. <none> 


265 






150.953 


ip0135 


tec 


pnUt!«p rimb irti \w & yjtcteit, 


Ononis) 


413" 






151,598 


!b0136 '..,u<\ 


nuUi-p i'n-u nalp ^in 


^(hone) 


199 






152,231 


ib0137 -vaJL 


i't"«!w^ firnbrt it pjol iin # 


Ipjone)^ 


202 




* 152.?54 


Tb0138 




paisdive nmbriaHtke protein 


*l<none) 


I 204 






155.428 !b01 39 


WrE 


probable outer membrane ponn protein t 


„l<none) 


866 






156.201 


Ib0140 


lecpD 


probable pUin thapsrene similar to PapD 


i (nmB) 


247 






! 150 233-tiPUI 


IvsW 


puiative fimbfisMIke protein 


< (none) 


195 






157,732 [S5T42 


ifOtK 


7,8-rjinVdro-6-h>^roKvmathyiptef}fi- pyrop. 


i{none) 


1 160 






159,093 ( &Q1 43 


\pQT)S 


poiyWj3olymerase I 


Unone) 


455 






160.112 Sb0l44 


tyacSB 


puiative. tRNA synthetase 


Unone) 


309 



















1 l_ocus f4ame 


Putative ID , | 


00114 


bvruvate dehydrogenase (decarboxyia^! 


••jb0115 


^rWate dehydrosehase (dlnydrbupoy../'- 






*ceEF ; . > pi 



ipyruvale derrydrogenase 



New 



Deietej 



^ fir bteins: 




^Abbreviation 




FIGURE 11 



WO 03/106998 



PCT/US03/18838 



12/18 



0i<5enelic Circuit Simufatioa - Madrfltevelopmenl 



ftle U:f Vjew Tools Help 




FIGURE 12 



WO 03/106998 



PCT/US03/18838 



13/18 



^{Genetic Circuit 




FIGURE 13 



WO 03/106998 



PCT/US03/18838 



14/18 



f&Genetie Circuit Simulation - Model Development 



Hie P'i'- View Tools Help ^. _ 



IS 



Model Content Browser 



jMou«5 |||N3ene fcrig* :j Protein index !^atS}on3n(4ek^ 

Content 



N«fc«oife 
Oaf* 



-Validalitm 



Mo<J*i 



Full Name 



| pyruvate del 



jb0114*b0i1S. 



., . -tem^&^'^^M',* 


! Locus Name | 


. PulatwelD ! 


&01 1 4 {pyruvate dehydrogenase t'decarboxyla,..' . 


bOjiS Jj^wSte deb^rogenase (dim'droii'po^..! 


IVo'teim 








pyruvate d\y°J^»6fL.„J 




< | ? Mew [| Delete : "| 


1 . < , 






Proteus •*;•;• ., : " ■)■ s >- ■ 


Abbreviation j V. 


Official Name > | 



r — 




(g) 



FIGURE 14 



WO 03/106998 



PCT/US03/18838 



15/18. 



g^Geneac Circuit Simulation - Model Development 



,ioixi 



FHe Ur- View Tools Htiip ■ \ \L^m I j?. ,. "'.V L\ 

f®SI®iSS^l i Qltoati M ortelj ^fttial; E«ti«tlchla ceil K.1J£.*Op«n EdHian* - 67RAVNIK ♦ { 1 1 -43 AM ) 



M»4tl 
Centtnf 



StfU«t»H%. 



Modtf 
ValldaHevn 



; Gene fewtex Protein Index 



Full Name 






gig g 








Locus Name V ->.'.•: 


Putative ID '! 


bdi 1 4 pyruvate d ehydroQenase (decarbc«j»)3„ >; : 


bans * pyruvate dehyto^ • 

P-'lftOttir ' . V : ,", ; j 


*eeEF ' ; ;A 




;pyruvate dehydrogenase ^ 




| New || 








j Abbreviation 1 ' ; 


Officialisms j 


I 




i „ IU.J . 












FIGURE 15 



WO 03/106998 



PCT/US03/18838 



16/18 



uit Stmulation - Modef Development 



F86 ViZW TOOJS HtiJp 



Content 



Mode) 
Piop*rtits' 



gj Content Browser ' Y-'r y V' ... 



: b 



fGenefrtdex,; Protein Index |:ReactiQnj«0px;j 



*ceEF. 



p ruvafe dl hydropenase bOH^ft^l'Vg " 



PYRDH 



Locus Name ; 



rate dehydrogenase {decatboxyte.,.1 



pyruvate dehydrogenase 



DfficiaiNarng 



*ceEF jpynrtate dehydrogenase! , 

pyrohT 



]{ , , •Pyruvate-dehydrogenase 
tea + : 'riad + pyr-> accoa + 602 + h + nadii • ;£ 




FIGURE 16 



WO 03/106998 



PCT/US03/18838 



17/18 



^Genetic Circuife Simulation - Model Development 



Wew Hetp 



iiiH 



l*gjLoadjMtodel[! .ftftttJtefc &sxM*mn* *>u w.i2.-<dp«« Edition* - sTR*VNiK-c<»i2W2oa2 oa^s:p^3 



.JS 

Content 

It 

StrueMJi* 



M#d*f 



Eg MQUetComeirt Browser ^ - LIU 



Gene Jmfax Pio+em Index if 



Ref! 


«*bre¥i..| 


!'4ame: 


Equation 


Confid.4 Subsystem '■ 


Region 


j 




Acn s 


Aeetalydehyde dehydrog,.. j 


accoa + <2) h + (2) nadh co... 


A 1 - 










Adenylate kinase 


amp + atp «~^£2) adg_ 


2 




;■ 




AKODH 


2-Osogluteraie dehydrog... 


2ogfutr * coa + nad -* co2 * h «v 


4 




| 




ATFASE 


Enoiase 


adgjcj «■ (3) hlel * pile} -> atplcj .,^ 


— It r 








CfTSYM 


Citrate synthase 


accoa + h2o * oaa -> cit + coa 


0! 




> 






Cytochrome oxidase 


<4> hid + olfcT + (855wnaR© , r . 


6i 








FADDH 


FAB dehydrogenase 


fadh ♦ quinone-> fad + quinon... 


ol - 




■:' 




FBP 


Fructose-bisphosphalase 


DHSbp ♦ h2o ~> Df6p + pi 


oi 








FRD 


Fumarate reductase 


fadh «-fum~> fad + sue 










FUMR 


Fumarase 


(urn + h2o <=> Lrnal 


of 1 






H I03PDH 




DjijyadSo + nad + pi > 3pfl1y~ 










ICDH 


Iso citrate dehydrogenase 


Isocit ♦ riadp-* 2og!utr + co2 * .J 


0 








isour 


isocrtrate lyase 


isocit->gryx+$uc 










MALS.. 


Malate synthase 


accoa + jjly> + h2o -» Ln ;d - to i 


0 








MDH 


Maiate dehydrogenase 


Lmal * nad <=-> h + nadh + oaa 


0; 








NADH... 


NADH dehydrogenase 


J3)hjc] + nadj^c} * quinoneicj-*..^ 


0! 








NADT,., 


WAD Transhvdrogenase 


nad + riatiph -> h + nadh nadp 


Of 




:; : 




pck 


Phosphoenoipyruvats . 


ato * oaa adp + co*2 *■ pip 


rSZIi 








PFK 


V.-r rWsphofructokin ice 
Phosphogiuconate dehy L 


DfBp + atp--> Oftebp + adp 
"epgicn ♦ nadp DrldtiSp + co2.„ 


0 








FODH 

PGK 


PhospboglyeerateWnase 


Spgiycp + adp <=> 3pg + atp 


0 








POM 
PPC 


FbosnhoofrceratemuW ; 3pg 2pg 


0? 








Phophoenoipyruyate can,. 


1 co2 + h2o + pep -> oaa * pi 


Oi 








PPS 


Phosphoenoipytuvaie sy... 


atp + h2o + pyr-^ amp + pep + pi 


F3fc== 




'■ ''>' 




PTA 


Ph< sphotraniaeetytase 


accoa + pi <=> actp + coa 










PYK 


Pyruvate kinase 


adp *j5ep-* atp^pyr 


01 








PYRDH 


iPjruvate dehydrogenase .. 


coa + nad + pyr ~> accoa ♦ co2 










RPE 


PhosphoriDufose eplmer.. 


Dribusp «sms» owisp 


Oi 




I 




RPl 


Phosphopentoisomerase 


Orihuop <=> OrBp 


1 0! 








S8UCOH 


C irrinr.1t eel sfdrogenase 


fad + sue ~> fadh + turn 


rz..lt 3" 








sues... 


Succinyi-coA synthetase 


* adp_f pi + succoa <~> atp ♦ co.. 


r oj 








TAL 


Transaidoiase 


Dah'ad3p * Ds?p «™> Dfi4p * ... 








I men 


iTransketaiase 1 


DrSp + Ox^lSp «ss> Dgryad3p + 


b| 






ITKT2 


fransketoiase 2 


De4u + D^'lon Dfop ♦ DQt/. 


6] 








ESlZ 


barbon dioxide transports 


dhap <™> Dglyad3p 


0 








TtranI 


{iC62te)<^> C021C] 


' o 






ITRAN... 


Ethanoi transporter 


etohje] etohlc) 


0 








[TRAM- 


Oxygen transporter . 


o2(ejj^> o2fc] 


0 








1TRAN... 


Phosphate transporter 


pije] *~~> pllcj 


1 0 








jTRAN^ 


succinaie Transporter 


suctel sucfc] 


j ° ! 






< 




• : ™ l> 





Beacitcm: ^ 

'agon" ~ , Aconttas©_ 



[Confidence Details t Notes 1 



\ 


Hot 












Hint | 


- am 





biochemical ,4 ' S { G r J> D 

feenetic Date ; .'B • .T G '3 

Sequence D,.. . 155 . .□ D O 

»hystologlca.,i B 1 O?""" 

Modeling Data SEi ; '□ □ ' 




FIGURE 17 



WO 03/106998 



PCT/US03/18838 



18/18 



Genetic Circuit Simulation - Model DeveJoprnent 



View Toots Help 



=asl*! 



gg Model Content Browser j-'.,;.. " r: :: ;- ' feg , -A v.; ; ^ ; : : • \ A , , AAy ±,A 



Content 



Exchange 
Reactions 

Strurtui* 

VaUdiriloa 



Motel V 



Gum Index PrtrtemJm.cx j Reaction Index j 





f ] AH 












Rett. 


^bbievj.j 


Name ' ) \ 


Equation 


Confid.. 


Subsystem ft 


, Region 


! 




ACT! — i Acetalydehyde dehydrog... 


accoa + 0 h + <2) nadh to... 


4 










ADK j 


Adenylate kinase 


amp + atp <=> (2) adp 












akodh; 


2-Oxosiuterate dehydrog... 


2oglutr + coa + nad «-» co2 + h +■,,. 


4. 






ATPASlfEnolase 


adp|c3 * <3) h[e] * pljcl-* atplcj ... 










OITSVW 


Citrate synthase 


accoa ♦ h2o * oaa -> cit * coa 


Oi 






ICYTOX 


Cytochrome oxidase _j 


(4) hfcj + o26^ + (2) quinoneRfcl ... 


0 








Ifaddh 


FAD dehydrogenase 


tedh-+ quinb'ne~* fed + qulnon... 


D 








1FBP 


Fmctase-bisphosphatasB lOftSbp * K2o ~> DP6p ♦ pi 


0 








1 FRD 


Fumarate reductase 


fadh ♦ fum -> fad + sue i 


0 








03PDH * 


Fu m ajase^ „ limw ,,. 


fUm* h2q<==* tnial 

Dar?ad¥p^*;ni¥+ pi 3pgife > 


0 
0 


tfnt< .. 1 -„ [ .„ t -.- ri - 


l\ 






ICDH ^ 


Isocitrate dehydrogenase 


Isocit * nadp ~* 2og)utr + co2 * ... 


0^ 








ISOLV 


Isocttrate lyase 


jsocu-> ct)yx*suc 


0* 








MALS 


Malate synthase 


accoa * glyx + h2o ~> Lmai + coa 


0' 








MDH 


Mai ate dehydrogenase 


tmal * nad h + nadh + oaa 


6 




K 




NADH... 


NADH dehydrogenase 


<3) h(c] * nadh{cj * ouinonelcj - 


Oi 








NADT . 


IMAD Transhydrogenase 


nad + nadph h + nadh + nadp 


0 










PCK 




atp + oaa ~> adp + co2 + pep' 


D 










PFK 


$-Phosphofruotok5nase 


t?!v,n - atp- r ifuh,n + adp 


0 










PGDH 


Phosphogluconate dahy... 


6pgicn + nadp ~* DrldUSp + eo2. c 


0 










POK 


Phosphogiycerateklnase 


3pgiycp + adp <=^=^ 3pg + atp 


OJ 








POM 




3prj <™> 2pg 


m 








PPO 


Fhlrphoerwipyiwatecar... 


co2 + TS'o* : pep ~> oaa + pi 


0 




.. 




PPS 


Phospnoenolpyruvate sy... 


atp + h2o ♦ pyr~» amp + pep + pi 


a 








PTA 


Phosphotransaeetylase 


accoa + Pi <sssi * actp + coa 


i 0 










iPYK 


Pyruvate kinase 


adp +j3ep~» atpjf- pyr 


0 






, 




PYROH 


jpjgj^tedehydrogenase 


coa + nad + pyr ~> accoa * co2 .„ 










;RPE 


Phosphortbulose epimir.. 


ormusp ox^isp 










RPl 


Phosphopentolsomerass 


|Dribu5p Drop 


- * 








SUCDH 


Succinate dehydrogenase 


ifid + sue -> fadh + turn 








sues... 


fsuccirwi-coA synthetase 


adp + pi + succoa alp + co.., 


6i 








TAL 


Transaldolase 


DalyadSp * 057 p Oe4p ♦ ... 


01 








TKT1 


Transketoiase 1 


DrSp + Oxyi6p <sas> Dfilyad3p + . 


01 








TtCT2 


^Transkelolase 2 


1 Oe4jj * Dwisp <sss5> D^p *• Dgty.. 


» 0 










TPI 


Triose-phosphate Isorne.. 


^dhap <—> Dglyad3p 


0 










TRAN,., 


carbon dioxide transporter co2le] <~~> co2|c] 


0 








TRAN.„ 


Ethanol transporter 


1 etoh|e} «"> elohlcj 


0^ 








TRAN 


Oxypen transporter . 


U2{eJ o2 (c) 


0 








'TRAlvT 


Vphosphate transporter 


r p|fe] p\[t\ 


0? 








TRAN... 


Succinate Transporter 


sucie] sucfej 


J01 ,j m 








■' 













offOH A v. f^sphaTede^rpsenase ] 

bgiyadap + nad * pf-^* 3pfilycp * h ^«^h 



Q»tnd»«c« ?o.o 



; Confidence Details |^M^s;;| 




ieneticData ; Pj ,Cr . p ; j Q \- O? 
Sequence D... , M t P ' \ Q □ ; P 

' |p 'Hp 



Physlologica..] p r P. 
Modeling Data[ I'D 









hf416 








FIGURE 18 



INTERNATIONAL SEARCH REPORT 



International application No. 
PCT/US03/18838 



A. CLASSIFICATION OF SUBJECT MATTER 

EPC(7) : G01N 33/48 
US CL : 702/19 
According to International Patent Classification (IPC) or to both national classification and IPC 

B. FIELDS SEARCHED 

Minimum documentation searched (classification system followed by classification symbols) 
U.S. : 702/19, 20; 703/2; 706/13 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 



Electronic data base consulted during the international search (name of data base and, where practicable, search terms used) 
Please See Continuation Sheet 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category * 


Citation of document, with indication, where appropriate, of the relevant passages 


Relevant to claim No, 


Y 


US 6,351,712 Bl (STOUGHTON et al.) 26 February 2002 (26.02.2002), see Detailed 


1, 7-17, 23-34, 42, 46- 




Description of the Invention. 


51, 57-58, 60, 69-76 


X 


US 6,132,969 A (STOUGHTON et al.) 17 October 2000 (17.10.2000), see Abstract and 


1-76 




Detailed Description. 




Y 


US 6,326,140 Bl (RINE et al.) 04 December 2001 (04.12.2001), see Detailed Description 


1-76 




of the Invention. 




Y 


US 6,329,139 Bl (NOVA et al.) 11 December 2001 (11.12.2001), see Summary of the 


1-76 




Invention. 




Y 


US 6,370,478 Bl (STOUGHTON et al.) 09 April 2002 (09.04.2002), see Detailed 


1-76 




Description of the Invention. 




X 


ROMERO et al. Nutrient-Related Analysis of Path way /Genome Databases. Pacific 


1-11, 17-28, 74-75 




Symposium on Biocomputing. 2001, pages 471-482, see Abstract, Introduction, and 




Y 


Algorithm. 


34-44, 51-65, 76 


X 


JUTY et al. Simultaneous modelling of metabolic, genetic and product-interaction networks. 


Ml, 17-28, 74-75 




Briefings in Bioinformatics. 2001, Volume 2, Number 3, pages 223-232, see Abstract and 




Y 


Information From the Web and Data Handling, and Creation of Model Files. 


34-44, 51-65, 76 



Further documents are listed in the continuation of Box C. 



□ 



See patent family annex. 



* Special categories of cited documents: 

"A" document defining the general state of the art which is not considered to be 
of particular relevance 

"E" earlier application or patent published on or after the international filing date 

"L" document which may throw doubts on priority claim(s) or which is cited to 
establish the publication date of another citation or other special reason (as 
specified) 

"O" document referring to an oral disclosure, use, exhibition or other means 

"P" document published prior to the international filing date but later than the 
priority date claimed 



later document published after the international filing date or priority 
date and not in conflict with the application but cited to understand the 
principle or theory underlying the invention 

document of particular relevance; the claimed invention cannot be 
considered novel or cannot be considered to involve an inventive step 
when the document is taken alone 

document of particular relevance; the claimed invention cannot be 
considered to involve an inventive step when the document is 
combined with one or more other such documents, such combination 
being obvious to a person skilled in the art 

document member of the same patent family 



Date of the actual completion of the international search 
22 September 2003 (22.09.2003) 


Date of mailing of the international search report 

04 NOV 2003 


Name and mailing address of the ISA/US 
Mail Stop PCT, Attn: ISA/US 

Commissioner for Patents / 
P.O. Box 1450 

Alexandria, Virginia 22313-1450 
Facsimile No. (703)305-3230 


Aumorizeja offlc^^? ^ S e . 
Telephone No. 703-308-0196 



Form PCT/ISA/210 (second sheet) (July 1998) 



INTERNATIONAL SEARCH REPORT 



International application No. 
PCT/US03/18838 



Box I Observations where certain claims were found unsearchable (Continuation of Item 1 of first sheet) 

This international report has not been established in respect of certain claims under Article 17(2)(a) for the following reasons: 

1. Q Claim Nos.: 

because they relate to subject matter not required to be searched by this Authority, namely: 



I | Claim Nos.: 

because they relate to parts of the international application that do not comply with the prescribed requirements to such 
an extent that no meaningful international search can be carried out, specifically: 



3 - EZ] Claim Nos.: 

because they are dependent claims and are not drafted in accordance with the second and third sentences of Rule 6.4(a). 

Box II Observations where unity of invention is lacking (Continuation of Item 2 of first sheet) 

This International Searching Authority found multiple inventions in this international application, as follows: 
Please See Continuation Sheet 



As all required additional search fees were timely paid by the applicant, this international search report covers all 
searchable claims. 

2. 1 1 As all searchable claims could be searched without effort justifying an additional fee, this Authority did not invite 

payment of any additional fee. 

3. 1 1 As only some of the required additional search fees were timely paid by the applicant, this international search report 

covers only those claims for which fees were paid, specifically claims Nos.: 



4* 1 1 No required additional search fees were timely paid by the applicant. Consequently, this international search report is 
restricted to the invention first mentioned in the claims; it is covered by claims Nos. : 

Remark on Protest [ [ The additional search fees were accompanied by the applicant's protest. 

I 1 No protest accompanied the payment of additional search fees. 



Form PCT/ISA/210 (continuation of first sheet(l)) (July 1998) 



INTERNATIONAL SEARCH REPORT 



PCT/US03/18838 



C. (Continuation) DOCUMENTS CONSIDERED TO BE RELEVANT 



Category ' 



X 
Y 



Citation of document, with indication, where appropriate, of the relevant passages 



MOSZER. The complete genome of Bacillus subtilis: from sequence annotation to data 
management and analysis. FBBS Letters. 1998, Volume 430, pages 28-36, see pages 28-35. 



JENSSEN et al. A literature network of human genes for high-throughput analysis of gene 
expression. Nature Genetics. 2001, Volume 28, pages 21-28, see Introduction. 

BECKERS et al. Large-scale mutational analysis for the annotation of the mouse genome. 
Current Opinion in Chemical Biology. 2001, Volume 6, pages 17-23, see Gene-driven 
aproaches and Phenotype-driven approaches. 

HARDISON et al. Globin Gene Server: A Prototype E-Mail Database Server Featuring 
Extensive Multiple Alignments and Data Compilation for Electronic Genetic Analysis. 
Genomics. 1994, Volume 21, pages 344-353, see Abstract. 



Relevant to claim No. 



1-11, 17-28, 74-75 



34-44, 51-65, 76 
1, 17, 34, 51, 74-76 

1-33, 74-75 



1-3, 17-19, 34-38, 51- 
53, 74-76 



Form PCT/ISA/210 (second sheet) (July 1998) 



INTERNATIONAL SEARCH REPORT 



PCT/US03/18838 



BOX II. OBSERVATIONS WHERE UNITY OF INVENTION IS LACKING 

This application contains the following inventions or groups of inventions which are not so linked as to form a single inventive 
concept under PCT Rule 13.1. In order for all inventions to be searched, the appropriate additional search fees must be paid. 

Group I, claims 1-33 and 51-75, drawn to a computer-implemented process and system for constructing a scalable 
output network model of a bioparticle. 

Group II, claims 34-50 and 76, drawn to a computer implemented process for self-optimizing a network model of a 
bioparticle. 

The inventions listed as Groups I-II do not relate to a single inventive concept under PCT Rule 13.1 because, under PCT Rule 
13.2, they lack the same or corresponding special technical features for the following reasons: The special technical feature of 
Group I is the production of a mathematical description of reactant fluxes. The special technical feature of Group II is the 
identification of an ameliorating network reaction component capable of augmenting the competence of the network model. 

Clearly, these two Groups with their respective technical features are distinct .from each other. Thus, in summary, the 
inventions of Groups I and II are not so linked under PCT Rule 13.1. 



Continuation of B. FIELDS SEARCHED Item 3: 

EMBASE, BIOSIS, SCISEARCH, MEDLINE, PUBMED, WEST searching the following terms: computer, network, process, model, 
bioparticle, compound, database, gene, open reading frames, reaction, transform, connectivity, mathematical, system, annotation, 
phenotype, self-optimizing, data structure, output, demand flux, stoichometric coefficients 



Form PCT/ISA/210 (second sheet) (July 1998) 



